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Abstract: Giant micelles are elongated, polymer-like objects created by the self-assembly 
of amphiphilic molecules (such as detergents) in solution. Giant micelles are typically flexible, 
and can become highly entangled even at modest concentrations. The resulting viscoelastic 
solutions show fascinating flow behaviour (rheology) which we address theoretically in this 
article at two levels. First, we summarise advances in understanding linear viscoelastic 
spectra and steady-state nonlinear flows, based on microscopic constitutive models that 
combine the physics of polymer entanglement with the reversible kinetics of self-assembly. 
Such models were first introduced two decades ago, and since then have been shown to 
explain robustly several distinctive features of the rheology in the strongly entangled regime, 
including extreme shear-thinning. We then turn to more complex rheological phenomena, 
particularly involving spatial heterogeneity, spontaneous oscillation, instability, and chaos. 
Recent understanding of these complex flows is based largely on grossly simplified models 
which capture in outline just a few pertinent microscopic features, such as coupling between 
stresses and other order parameters such as concentration. The role of 'structural memory' 
(the dependence of structural parameters such as the micellar length distribution on the flow 
history) in explaining these highly nonlinear phenomena is addressed. Structural memory 
also plays an intriguing role in the little-understood shear-thickening regime, which occurs 
in a concentration regime close to but below the onset of strong entanglement, and which is 
marked by a shear-induced transformation from an inviscid to a gelatinous stately 
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CHAPTER 1. RHEOLOGY OF GIANT MICELLES 



1.1 Introduction 



An amphiphilic molecule is one that combines a water-loving (hydrophilic) part or 'head 
group' with a with a water-hating (hydrophobic) part or 'tail'. The head-group can be 
ionic, so that the molecule becomes charged by dissociation in aqueous solution; or nonionic 
(but highly polar, favouring a water environment), in which case the amphiphile remains 
uncharged. Zwitterionic head-groups, with two charges of opposite sign, are also common. 
The hydrophobic tail is almost always a short hydrocarbon (though fluorocarbons can also 
be used); in some cases (such as biological lipids) there are two tails. The most important 
property of amphiphilic molecules, from the viewpoint of theoretical physics at least, is their 
tendency to self-assemble by aggregating reversibly into larger objects. The simplest of these 
is a spherical aggregate called a 'micelle' which in water h as t he hydrophobic tails sequestered 
at the centre, coated by a layer of headgroups; see Fig jl.li (In a nonaqueous solvent, the 
structure can be inverted to create a 'reverse micelle'.) 

For geometrical reasons, a spherical micelle is self-limiting in size: unless the amphiphilic 
solution contains a third molecular component (an oil) that can fill any hole in the middle, 
the radius of a micelle cannot be more than about twice the length of the amphiphile. To 
avoid exposing tails to water, it also cannot be much less than this; the resulting 'quorum' of 
a few tens of molecules for creation of a stable micelle leads to a sharp minimum in free energy 
as a function of aggregation number. This collective aspect to micelle formation causes the 
transition from a molecularly dispersed solution to one of micelles to be rather sharp; micelles 
proliferate abruptly when the concentration is raised above the 'critical micelle concentration' 
or CMC [1]. 

In a spherical micelle the volume ratio of head- and tail-rich regions is also essentially 
fixed: such micelles are favoured by amphiphiles with relatively large size ratios between head 
and tail. Suppose this size ratio is smoothly decreased, for instance by adding salt to an 
ionic micellar solution (effectively reducing the head-group size by screening the coulombic 
repulsions). The most stable local packing then evolves from the spherical micelle towards a 
cylinder; Fig jl.li (Proceeding further, it becomes a fiat bilayer; systems in which this hap- 
pens are not addressed here.) Allowing for entropy, the transition from spheres to cylinders 
is not sudden, but proceeds via short cylindrical micelles with hemispherical end-caps. The 
bodies of these cylinders smoothly increase in length as the packing energy of the body falls 
relative to the caps; the micelles eventually become extremely long. The law of mass action, 
which favours larger aggregates, means that in suitable systems the same sequence can be 
observed by increasing concentration at fixed head/tail size ratio (fixed ionic strength). 

Since the organization of amphiphiles within the cylindrical body is (in most cases) fluid- 
like, the resulting 'giant micelles' soon exceed the so-called persistence length, at which ther- 
mal motion overcomes the local rigity and the micelle resembles a flexible polymer chain. 
This crossover may or may not precede the 'overlap threshold' at which the volume occu- 
pied by a micelle - the smallest sphere that contains it - overlaps with other such volumes. 
Beyond this threshold, the chainlike objects soon become entangled but (unless extremely 
stiff) remain in an isotropic phase with no long-range orientational order. At very high con- 
centrations, such ordering can arise, as can positional order, giving for example a hexagonal 
array of near-infinite straight cylinders. Branching of micelles can also be important in both 
isotropic and ordered phases; a schematic phase diagram, applicable to many but not all 
systems containing giant micelles, is shown in Fig jl.li 



This article addresses primarily the isotropic phase of giant micelles, in a concentration 
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Figure 1.1: Upper panel: Schematic view of aqueous surfactant self-assembly. Lower panel: 
Schematic phase diagram for self assembly of ionic amphiphiles into giant micelles and related 
structures. The vertical axis represents volume fraction $ of amphiphile; the horizontal is 
the ratio Cs/C of added salt to amphiphile concentrations. Figure in lower panel from F. 
Lequeux and S. J. Candau, in Ref. [2], reprinted with permission. 
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range from somewhat below, to well above, the overlap threshold. This is the region where 
viscoelastic (but isotropic) solutions are observed. This regime merits detailed attention for 
two reasons. The first is that micellar viscoelasticity forms the basis of many applications, 
ranging from personal care products (shampoos) to specialist drilling fluids for oil recov- 
ery [2]. The second, is that as emphasised first by Rehage and Hoffmann [3], viscoelastic 
micelles provide a uniquely convenient laboratory for the study of generic issues in nonlin- 
ear flow behaviour. This is partly because, unlike polymer solutions (which they otherwise 
resemble), the self-assembling character of micellar solutions causes them to self-repair after 
even the most violently nonlinear experiment. (In contrast, strong shearing of conventional 
polymers causes permanent degradation of the chains.) Our focus throughout this review is 
on rheology, which is the science of flow behaviour. Although we will refer in many places 
to experimental data, we make no attempt at a comprehensive survey of the experimental 
side of the subject, nor do we describe applications areas in any detail. For an up-to-date 
overview of both these topics, the reader is referred to a recent book [4] , of which a shortened 
version of this article forms one Chapter. 

We s hall address the theoretical rheology of giant micelles at two levels. The first (in 
Section 11.41) is microscopic modelling, in which one seeks a mechanistic understanding of 
rheological behaviour in terms of the explicit dynamics — primarily entanglement and re- 
versible self-assembly — of the giant micelles themselves. This yields so-called 'constititive 
equations' which relate the stress in a material to its deformation history. Solution of these 
equations for simple experimental flow protocols presents major insights into the fascinat- 
ing flow properties of viscoelastic surfactant solutions, including near-Maxwellian behaviour 
(exponential relaxation) in the linear regime, and drastic shear-thinning at higher stresses. 
These successes mainly concern the strongly entangled region where the micellar solution is 
viscoelastic at rest; in this regime, strong shear-thinning is usually seen. There are however 
equally strange phenomena occurring at lower concentrations where the quiescent solution is 
almost inviscid, b ut beco mes highly viscoelastic after a period of shearing. These will also be 
discussed (Section 11.4.61) althought they remain, for the present, much less well understood. 

Microscopic models of giant micelles under flow generally treat the micelles as structure- 
less, flexible, polymer-like objects, albeit (crucially!) ones whose individual identities are 
not sustained indefinitely over time. This neglect of chemical detail follows a very successful 
precedent set in the field of polymer dynamics [5,6]. There, models that contain only four 
static parameters (persistence length, an excluded volume parameter, the concentration of 
chains, and the degree of polymerization or chain length) and two more dynamic ones (a 
friction constant or solvent viscosity, and the so-called 'tube diameter') can explain almost 
all the observed features of polymeric flows. Indeed, microscopic models of polymer rhe- 
ology arguably represent one of the major intellectual triumphs of 20th century statistical 
physics [7]. 

However, at least when extended to micelles, these microscopic constitutive models remain 
too complicated to solve in general flows, particularly when flow instabilities are present. 
(Such instabilities are sometimes seen in conventional polymer solutions, but appear far 
more prevalent in micellar systems.) Moreover, they omit a lot of the important physics, 
particularly couplings to orientational fields and concentration fluctuations, relevant to these 
instabilities. Therefore we also describe in Section 11.51 some purely macroscopic constitutive 
models, whose inspiration stems from the microscopic ones but which can go much further in 
addressing the complex nonlinear flow phenomena seen in giant micelles. These phenomena 
include for example "rheochaos", in which a steady shear deformation gives chaotically 
varying stress or vice versa. Our discussion of macroscopic modelling will take us to the 
edge of current understanding of these exotic rheological phenomena. 
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Prior to discussing rheology, we give in Section 11.21 a brief survey of the equilibrium 
statistical mechanics of micellar self-assembly. More detailed discussions of many of the 
static equilibrium properties of micelles can be found in [4]; we focus only on those aspects 
needed for the subsequent discussion of rheology. Another key component for rheological 
modelling is the kinetics of micellar 'reactions' whereby micelles fragment and/or recombine. 
These reactions are of course already present in the absence of flow, and represent the kinetic 
pathway whereby equilibrium (for quantities such as the micellar chain-length distribution) 
is actually reached. We review their properties also in Section [L2l 

In developing the equilibrium statistical mechanics and kinetic theory for giant micelles 
(Section II. 2p . we should keep in mind both the successes and limitations of the rheological 
theories that come later. Such theories, since their first proposal by one of us in 1987 [8] 
have had considerable success in predicting the basic features of linear viscoelastic relaxation 
spectra observed in experiments, and in inter-relating these, for any particular chosen sys- 
tem, with nonlinear behaviour such as the steady-state dependence of stress on strain rate. 
These dynamical models take as input the micellar size distribution, stiffness (or persistence 
length) and the rate constants for various kinetic processes that cause changes in micellar 
length and topology. Such inputs are theoretically well defined, but harder to measure in 
experiment. Nonetheless, there are a number of 'primary' predictions (such as the shape 
of the relexation spectrum, and the inter-relation of linear and nonlinear rheological func- 
tions; see Section n. 4. 3p for which the unknown parameters can either be fully quantified, or 
else eliminated. As an aid to experimental comparison, it is of course useful to ask how the 
rheological properties should depend on thermodynamic variables s uch as surfactant concen- 
tration, temperature, and salt-levels in the micellar system (Section ll.4.4p . But in addressing 
these 'secondary' issues, the dynamical models can only be as good as our understanding of 
how those thermodynamic variables control the equilibrium micellar size distribution, per- 
sistence length, and rate constants, as inputs to the dynamical theory. In many cases this 
understanding is only qualitative, so that these 'secondary' experimental tests should not be 
taken as definitive evidence for or against the basic model. 



1.2 Statistical Mechanics of Micelles in Equilibrium 



In line with the above remarks, we focus mainly on those aspects of equilibrium self-assembly 
that can affect primary rheological predictions. Most of th e thermodyn amic modelling can 
be addressed within mean-field-theory approaches (Sections ll.2.1f[I.2.3l) . although more ad- 
vanced treatments show various subtleties that still await experimental clarification (Section 
ll.2.4p . In Section [1.2.51 we turn to the kinetic question of how micelles exchange material 
with one another within the thermal equilibrium state. 



1.2.1 Mean Field Theory: Living Polymers 

In typical giant micellar systems the critical micelle concentration (CMC) is low - of order 
10~^ molar for CTAB/KBr, for example. (CTAB, cetyltrimethylammonium bromide, is a 
widely studied amphiphile. In what follows, we do not expand the acronyms for this or other 
such materials as their chemical formulas are rarely of interest in our context. KBr is, as 
usual, potassium bromide, added to alter the head-group interactions.) As concentration is 



6 



CHAPTER 1. RHEOLOGY OF GIANT MICELLES 



raised above the CMC, uniaxial elogation occurs and soon micelles become longer than their 
persistence length Ip. This is the length over which appreciable bending occurs [5]; once 
longer than this, micelles resemble flexible polymers. Persistence lengths of order 10 - 20 nm 
are commonplace, though much larger values are possible in highly charged micelles at low 
ionic strength. 

As concentration is increased, there is an onset of viscoelastic behaviour at a volume 
fraction C usually identified with a n 'ove rlap' concentration C* for the polymers. (For 
problems with this identification, see 11.2.31 below.) Above C*, the wormlike micelles are in 
the so-called 'semidilute' range of concentrations [5] - overlapped and entangled at large 
distances, but well separated from one another at scales below ^, the correlation length or 
'mesh size'. In ordinary polymer solutions in good solvents, the behaviour at scales less than 
^ is not mean-field-like but described by a scaling theory with anomalous exponents [5] . We 
return to this in Section II. 2. 4^ but note that these scaling corrections become small when 
the persistence length of a micellar cylinder is much larger than its diameter, giving modest 
values for a dimensionless 'excluded volume parameter' w [5,6]. Therefore, a mean-field 
approach - in which excluded volume interactions are averaged across the whole system 
rather than treated locally - captures the main phenomena of interest, particularly in the 
regime of strong viscoelasticity at C > C*. 

The simplest mean field theor y [9, 10] assu mes that no branch-points and no closed rings 
are present (rectified in Sections 11.2.21 [1.2. 3p . and ascribes a free energy E /2 to each hemi- 
spherical endcap of a micelle relative to the free energy of the same amount of amphiphilic 
material residing in the cylindrical body. Denoting by c{N) the number density of aggregates 
containing N amphiphiles or 'monomers', the mean field free energy density obeys 

PF = J2c{N)[\nc{N)+pE] + Fo(0) (1.1) 

N 

Here (3 = l/ksT; the term in E counts two end-caps per chain, and the cine piece comes 
from the entropy of mixing of micelles of different lengths. Within a mean-field calcula- 
tion, these are the only terms sensitive to the size distribution c{N) of the micelles; the 
free energy (including configurational entropy) of the cylindrical sections, alongside their 
excluded- volume interactions and all solvent terms, give the additive piece Fo{(f)) which de- 
pends only on total volume fraction (p. (It may also depend on ionic strength and related 
factors.) The volume fraction obeys 

(P = VoC = VoY.Nc{N) (1.2) 

N 

where vq is the molecular volume of the amphiphiles and C their total concentration. 
Minimizing (11. ip at fixed gives an exponential size distribution 

c(A^) oc exp[-A^/iV] ; N ^ (f)^/"^ exp[(3E/2] (1.3) 

The exponential form in each case is a robust result of mean field theory. The ^-dependence 
in the second equation is also robust (it follows from mass action), but can be treated 
separately from the much stronger exponential factor only so long as parameters like E 
and Vq are themselves independent of concentration. (In ionic systems this is a strong and 
questionable assumption.) The formula for N as written in (11.31) suppresses prefactoral 
dependences on VQ, lp and Oq, where Oq is the cross-sectional area of the micellar cylinders; 
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these are absorbed into our definition of E. So long as oq is constant, tlien exactly the same 
functional forms as in fll.Sp control c{L) and L, where L oc is the contour length of a 
micelle. Within mean field, L in turn controls the typical geometric size R (usually chosen 
as either the end-to-end distance, or the radius of gyration) of a micelle via R^ ~ Lip. This 
is the well-known result for gaussian, random- walk chain configurations [5]. 

We can now work out, within our mean- field approach, the overlap concentration C*, 
or overlap volume fraction 0* = C*vo. For a micelle of the typical contour length L we 
have R ~ n^/'^lp where n = L/lp is the number of persistence length it contains; this obeys 
nZpOo/fo = A^. The total volume of amphiphile within the region spanned by this micelle is 
Nvq and the volume fraction within it therefore (p ^ Nvq/R^. At the threshold of overlap, 
this (f) equates to the true value 0*; then eliminating via (11.31) gives 



For typical cases the dimensionless pre-exponential factor is smaller than unity, but nonethe- 
less a fairly large E is required if (p* is to be below, say 1%. The regime of long, entangled 
micelles usually entails scission energies E of around 10 — 20A;bT; in practice, experimental 
estimates of (p* (best determined by light scattering) are often in the range 0.05-5% [11]. The 
scission energy E of course depends on the detailed chemistry of the surfactant molecules 
and this (alongside micellar stiffness or persistence length) is one of the main points at which 
such details enter the theory. Very crudely, one can argue that doubling the mean curvature 
of a micellar cylinder to create an end-cap must cost about ksT per molecule in the end- 
cap region. (If packing energies were much higher than this, one would expect a crystalline 
rather than fluid packing on the cylinder, which is not typically observed, at least at at room 
temperature.) This gives E ~ nksT where n is the number of molecules in two endcaps. 
Within a factor two, this broadly concurs with the range 10 — 20/cbT stated above. More 
precise theoretical estimates also concur with this range, although values well outside of it 
are also possible for atypical molecular geometries, e.g. fluorosurfactants [12]. 

Th e region around (p* is where spectacular shear-thickening rheology occurs (see Section 
ll.4.6p . In ionic micellar systems without excess of salt, the strong dependence of lp,E and 
other parameters on cp itself i n th is region means that the simple calculations leading to 
(11. 3p . and hence the estimate (II. 4p . are at their least reliable. More detailed theories, which 
treat electrostatic interactions explicitly, give a far stronger dependence of L on and also 
a narrower size distribution for the micelles [13]. The overlap threshold (p* itself moves to 
higher concentration due to the electrostatic tendency to stabilise short micelles. 



1.2.2 Role of Branching: Living Networks 

The above assumes no branching of micelles. A mean-fleld theory can in principle be formu- 
lated to deal with self-assembled micellar networks having arbitrary free energies for both 
end caps and branch points [14]. This is, however, somewhat intractable for the general 
case. Fortunately things simplify considerably in the branching-dominated limit; that is, 
when there are many branch-poi nts p er end-cap. For branching via z-fold 'crosslinks' (each 
of energy Eb) one has, replacing (II. ip . the following mean-fleld result [14]: 



(1.4) 



(3F = Y^ c{N) [In c{N) + (3Ei,z-^] + 2{z-^ - l)C ln(2C) + Fo(0) 



(1.5) 



N 
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where C is as defined in (11 ■2p . and c{N) is now the concentration of network strands con- 
taining amphiphiles. To understand this resuh, note that the first logarithmic term is 
the translational entropy of a set of disconnected network strands. The second such term 
estimates the entropy loss on gathering the ends of these strands locally to form z-fold junc- 
tion points. T he te rm in Eh counts the energy of these junctions and -Fo(0) has the same 
meaning as in (11. ip . The value of z most relevant to micelles is z = 3, since for a system 
whose optimal local packing is a cylinder, a three-fold junction costs less in packing energy 
than z > 3. Low z is also favoured entropically: to create a four- fold junction one must fuse 
two three- fold ones with consequent loss of translational entropy along the network [14]. 

Minimizing (II. 5p to find the equilibrium strand length distribution, one finds this again to 
be exponential, with mean strand length L ~ exp[/9£'b]. This result applies whenever 

the geometric distance between crosslinks, A ~ {LlpY^"^ greatly exceeds the geometric mesh 
size ^, which within mean field theory obeys ^ ~ ag/f/'/p. This situation of A ^ ^ is 
called an 'unsaturated network' [14] and arises at high enough concentrations (0 ^ 0'^'^* ~ 
VQ'^ex'p[—PE/(3 — z/2)]). For < 0'"'* one has a 'saturated network' with A ~ ^. At low 
enough this saturated network can show a miscibility gap, where excess solvent is expelled 
from the system rather than increase ^ which would sacrifice branch point entropy [14]. 

The rheology of living networks (see Section 11.4.51) should differ strongly from the un- 
branched micellar case. Such a regime has been identified in several systems, primarily 
cationic surfactants at relatively high ionic strength [15-17]. These accord with the expected 
trend for curvature packing energies: adding salt in these systems stabilizes negatively curved 
branch-points relative to positively curved end-caps [16]. 



1.2.3 Role of Loop Closure: Living Rings 

We have assumed in (ll.ip that rings do not arise. A priori, however, there is nothing to stop 
micelles forming closed rings. Moreover, for unbreakable polymers (at least) the effect on 
rheology of closing chains to form rings is thought to be quite drastic [18], so this assumption 
merits detailed scrutiny. It turns out to be satisfactory only when E is not too large, so that 
0* in (II. 4p lies well above a certain volume fraction (p^""^, defined below, which signifies a 
maximal role for ring-like micelles. 

From (ll.4p . as E oo, the overlap threshold 0* for open micelles tends to zero. In 
this limit, there is formally just a single micelle, of macroscopic length. This corresponds 
to an untenable sacrifice of translational entropy which is easily regained by ring formation. 
To study this, let us set ^ oo so that no open chains remain, but allow rings with 
concentration Cr{N). Then, to replace (II. ip . one has [19] 

f3F = J2cr{N)[lnCr{N) + (3fr{N)] + Fo(0) (1.6) 

N 

where fr{N) = —kBTln{Zr), and Z,. is the configurational free energy cost of ring closure. 
Put differently, E — fr{N) is the total free energy cost of hypothetically opening a ring, 
creating two new endcaps but gaining an entropy A;^ In . The latter stems both from the 
number of places such a cut could occur, and the extra configurations made available by 
allowing the chain ends to move apart. 

For gaussian (mean-field-like) chains in three dimensions, it is easly shown that Z^. = 
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\]\[ 5/2 j-gj^ where A is a dimensionless combination (as yet unknown [11]) of ao,fo,^p- Min- 
imizing (11.61) at fixed 4> tlien gives 

Cr{N) = AAT-^/^e'^^ (1.7) 

wliere ft is a. cliemical-potential like quantity. Interestingly, this size distribution for rings 
shows a condensation transition. That is, for /i > the volume fraction 0^ = VoJ^N^'^riN) 
is divergent, whereas for /i < it can apparently be no greater than 

oo 

N=Nmin 

This limiting value of 0^ depends not only on A but on Nmin, which denotes the smallest 
number of amphiphiles that can make a ring-shaped micelle without prohibitive bending 
cost. (Such a micelle must presumably have contour length of a few times Ip.) 

This mathematical situation, in which there is an apparently unphysical 'ceiling' 0™"^ < 1 
on the total volume fraction of rings a system can contain, is reminiscent of Bose condensation 
[20]. It represents the following physical picture, valid in the E ^ oo limit. For < 0^^^^ 
one has the power law distribution of ring sizes in (11.71) . cut off at large N by an exponential 
multiplier (resulting from small negative fi). For > 05!*"^, one has a pure power law 
distribution of rings, in which total volume fraction 0™"^ resides; plus an excess volume 
fraction — 0^"^ which exists as a single 'giant ring' of macroscopic length. This giant ring 
is called the condensate; its sudden formation at = 0^^*°^^ represents a true phase transition. 

Obviously, an infinite ring is possible only because the limit E oo was taken; for any 
finite E, all rings with > as defined in (II. 3p . including the condensate, will break up 

into pieces (roughly of size A^). Indeed, if E is finite one can, within mean-field theory, 
simply add the chain and ring free energy contributions as 

PF = Y.c{N)[\nc{N)+f3E]+Y.Cr{N)[lnCr{N)+f3fr{N)]+Fo{(l)) (1.9) 

N N 

From this one can prove that the condensation transition is smoothed out for any E < oo 
[21]. Nonetheless, if E is la rge enough that the overlap threshold 4>*{E) obeying (11.41) 
falls below 0™"^ obeying (II. 8p . then the condensation transition of rings, though somewhat 
rounded, should still have detectable experimental consequences. These should mainly affect 
a (roughly) factor-two window in concentration either side of 0™"^. For -C 0™"^ there are 
no very large rings and hence limited opportunities for viscoelasticity. For ^ (pmax ^^iq 
volume fraction — 0^ of long chains exceeds that of rings, and the chains dominate. 

Because of uncertainty over the values of A and A^mm in (H-Sp and how these might depend 
on chain stiffness, ionic strength, etc., (f)^°'^ is one of the least well- char act arized of all static 
quantities for giant micelles. In fact, there is relatively little (but some [22]) experimental 
evidence for a ring-dominated regime in any micellar system, suggesting perhaps that, for 
reasons as yet unclear, 0^"^ li es wel l below the range of accessed in most experiments. 
However, as outlined in Section fl .4.61 below, a ring-dominated regime might explain some of 
the strangest of all the rheological data in the shear-thickening regime just below C [23]. 

Note that in an earlier review (Ref. [11]) an impression was perhaps given that ring- 
formation matters only within scaling theories (discussed next) but not at the mean-field 
level. This is true only if 05!*°^ is indeed small; in that case rings will only matter for very 
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large E, and micelles are of sufficient size that excluded volume effects, even if locally weak, 
are likely to give scaling corrections to mean-field. However, rings are not purely a scaling 
phenomenon: even in a strictly mean-field picture, 0^"^ is a well-defined quantity at which 
a ring-condensation phase transition arises in the limit of large E. 



1.2.4 Beyond Mean Field Theory 

As with conventional polymers [5], micelles can exhibit non-gaussian statistics, induced by 
excluded-volume interactions arising from the inability of two different sections of the micelle 
to occupy the same spatial position. In principle, this gives scaling corrections to all the 
preceding mean-field results, altering the various power law exponents that appear in equa- 
tions such as (11.3111.71) . As mentioned previously, however, micelles often have a persistence 
length Ip large compared to their cross-section and therefore tend to have relatively weak 
excluded volume interactions. Therefore, giant micellar systems can often be expected to lie 
in a messy crossover region between mean-field and the scaling theory. For completeless we 
outline the scaling results here (see [11,14,24] for more details), but without attempting to 
track dependences on parameters like /p,ao,fo- 

First, consider a system with no branches or rings. In such a system, the excluded volume 
exponent u ~ 0.588 governs the non-gaussian behaviour of a self-avoiding chain; R L^ [5]. 

This gives where d = 3 (the dimension of space). This leads to a scaling 

of the transient elastic modulus Go (defined in Section [L^ below) : (3Go ~ ~ 0^'^, which 
is proportional to the osmotic pressure 11 [5]. This differs from the simplest mean-field-type 
estimate which has Go ~ n ~ 0^. Second, one finds in place of (II. 3p 

c(L) ocexp[-L/I] ; L ~ 0^ exp[/3E/2] (1.10) 

where y = [1 + (7 — l)/{vd — l)]/2 ^ 0.6; here 7 ~ 1.2 is another standard polymer 
exponent [5]. In practice this is rather a modest shift from the result in ( II. 3p . 

In the presence of branch points, the important case remains 2 = 3. Here the mean field 
result for the mean network strand length L ~ 0~-^/^ exp[/3i?] becomes L ~ exp[/5ii^] 
with an exponent A ~ 0.74; the expression for A in terms of standard polymer exponents 
is given in Ref. [14]. Similarly the mean- field result 0'^'^* ~ g-2/3£;/3 becomes 0'^"* ~ ^-l^E/y 
with y ~ 0.56 [14]. Note that there is still the possibility of phase separation between a 
saturated network and excess solvent, even under good solvent conditions - and there is 
some experimental evidence for exactly this phenomenon [17]. 

Scaling theories in the presence of rings become even more complicated [19]; even the 
exponent z/, governing the local chain geometry, is slightly different in the ring dominated 
regime [24]. Near 0jy*"^ (which has the same meaning as before, but no longer obeys (11.81) ) 
there is a power-law cascade of rings, controlled by a distribution similar to (II. 7p . but with a 
somewhat larger exponent (2.74 instead of 5/2) [24]. This cascade of rings creates 'power law 
screening' of excluded volume interactions [25], causing u to shift very slightly downward [19]. 
Perhaps fortunately in view of these complications, the ring-dominated regime, if it exists 
at all, is poorly enough quantified experimentally that comparison with mean field theory is 
all that can be attempted at present. 
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Figure 1.2: The three main types of micellar reaction: Top, reversible scission; middle, 
interchange; bottom, bond interchange 



end 



1.2.5 Reaction Kinetics in Equilibrium 



Alongside the micellar length distributions addressed above, a key ingredient into rheological 
modelling is the presence of reversible aggregation and disaggregation processes (which we 
shall call micellar 'reactions'), allowing micelles to exchange material. We will treat these 
reactions at the mean-field level, in which micelles are 'well-mixed' at all times (so there 
is no correlation between one reaction and the next). Our excuse for this simplification 
is that, although deviations from mean- field theory are undoubtedly important in some 
circumstances and have been worked out in detail theoretically [26], there is so far rather little 

evidence that this matters in the strongly entangled regime (C > C) where viscoelasticity 
is primarily seen. And, although in the shear thickening region (C ~ C) it is quite possible 
that non- mean- field kinetic effects become important, there is so much else that we do not 
understand about this regime (Section [1.4.61) that a detailed discussion of correlated reaction 
effects would appear premature. 

We neglect branching and ring-formation in the first instance, and also distinguish reac- 
tions that change the aggregation number of a particular micelle iV by a small increment, 
AN ~ 1, from those which create changes AA^ of order A^ itself. The former reactions can of 
course lead to significant changes in micellar size over time, but as A^ increases, the timescale 
required for this gets longer and longer [27]. Unless the reaction rates for all reactions of the 
second (AA^ ~ A^) type are extremely slow, these latter will dominate for large aggregates. 
From now on, we consider only reactions with AA^ ~ A^, of which there are thre e ba sic types: 
reversible scission, end-interchange, and bond-interchange, as shown in Figure 11.21 

In reversible scission, a chain of length L breaks spontaneously into two fragments of size 
L' and L" = L — L' . (Note that the conservation law A^" + A^' = A^ really applies to A^ 
and not L as written here; but if we ignore the minor corrections to L represented by the 
presence of end-caps, the sum of micellar lengths is also conserved.) In thermal equilibrium 
the reverse process (end-to-end fusion) happens with exactly equal frequency; this follows 
from the principle of detailed balance [28]. If, for simplicity, we assume that the fusion rate 
of chains of lengths L and L' is directly proportional to the product of the ir co ncentrations, 
then the fact that detailed balance holds for the equilibrium distribution (11. 3p can be used 
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to reduce the full kinetic equations (as detailed, e.g., in Ref. [8]) to a single rate constant 
krs- This is the rate of scission per unit length of micelle, and is independent of both the 
position within, and the length of, the micelle involved [8]. More relevant physically is 

Trs = (KsL)-^ (1.11) 

which is the time taken for a chain of the mean length to break into two pieces by a re- 
versible scission process. Note that, by detailed balance, the lifetime of a chain-end before 
recombination is also Trs [8]. Moreover, solution of the full mean-field kinetic equations [29] 
shows that if L is perturbed from equilibrium, for example by a small jump in temperature 
(T-jump), L relaxes monoexponentially to equilibrium with a decay time Trs/2. (In fact this 
applies not only to L but to the entire perturbation to the micellar size distribution, which 
for this form of disturbance is an eigenmode of the kinetic equations [29].) The response to 
a nonlinear, large amplitude jump is also calculable [30]. These results allow Trs to be esti- 
mated fro m T -jump data [31], providing an important constraint on the rheological models 
of Section 11.41 

Turning to end-interchange, this is the process wh ere a 'reactive' chain-end bites into 
another micelle, carrying away part of it (Figure II. 2p . Assuming all ends to be equally 
reactive, and applying detailed balance, one finds that all points on all micelles are equally 
likely to be attacked in this way. There is, once again, a single rate constant kei, but now 
the lifetime of any individual chain end is l/kei4>, since the availability of places to bite into 
is proportional to 0. The lifetime of a micelle of the average length, before it is involved in 
an end-interchange reaction of some sort, is [29] 

Te, = (4A;e.0)-^ (1.12) 

In contrast to the reversible scission case, analysis of the full mean-field kinetic equations [29] 
shows that end-interchange is invisible in T-jump: for the specific form of perturbation that 
arises there, no relaxation whatever occurs by this mechanism. Beyond mean-field kinetics 
this would no longer hold, but there remains an important limitation to end-interchange in 
bringing the system to equilibrium. Specifically, end-interchange conserves the total number 
of micelles. Accordingly if a disturbance, whether rheological or thermal in origin, is applied 
that perturbs the total chain number J2n c{N) away from equilibrium, this will not fully 
relax until the time-scale Trs is attained, even if this is much larger than Tei [29] . In the mean 
time, the end-interchange process relaxes the size distribution c{N) towards the exponential 
form c{N) oc exp[—N/N] of (11. 3p : but with a nonequilibrium value of A^. This separation 
of time scales m ay lie at the origin of strange 'structural memory' effects seen in certain 
systems (Section 11.4.81 below) [23]. 

Note that since in our simple models the micellar energy is fixed by the number of end caps, 
conservation of micellar number in end-interchange reactions (and also bond-interchange, 
below) is tantamount to conservation of the total energy stored in such end caps. An 
energy-conserving processes cannot, unaided, relax a system after a jump in temperature. 
However, since E is really a free energy and the dynamics is not microcanonical, conservation 
of micellar number is perhaps the more fundamental concept in distinguishing interchange 
from reversible scission kinetics; in subsequent discussions, we take this view. 

Finally we turn to the bond-interchange process [32] in which micelles transiently fuse 
to f orm a four-fold link before splitting again into differently connected components (Figure 
II. 2p . This process, like end-interchange, conserves chain number. Indeed it does not even 
alter the identity of chain ends. Since, in entangled polymeric systems, stress relaxation oc- 
curs primarily at the chain ends, bond-interchange is far less effective than reversible scission 
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or end-interchange in speeding up the disentanglement of micelles (see Section [L^ . In fact, 
although a breaking time Tbi = { kbiL(f)) ~^ can be defined, this enters the rheological models 
differently from Trs or (Section II .4.21 below) . Bond interchange also allows chains to effec- 
tively pass through one another by decay of the four-fold intermediate, creating a somewhat 
different relaxation channel for chai n dise ntanglement and stress relaxation [33]. However 
(as previously discussed in Section 11.2.21) a transient four-fold link is likely to dissociate 
rapidly into two three-fold links. Such three-fold links are in turn the transition states of the 
end-interchange process. If these links disconnect rapidly, then the end-interchange process 
(which their decay represents) is probably dominant over bond interchange. If they do not 
decay rapidly, then it is likely that their existence cannot be ignored for static purposes; one 
has a branched system in equilibrium (see Section [1.2.21) . 

The reaction kinetics in branched micellar networks is far from easy to cast in terms of 
simple mean-field reaction equations, as studied, e.g., in Ref. [29] for unbranched chains. 
However, within such a network, alongside any bond- interchange reactions that are present, 
structural relaxation can still occur by reversible scission or end-interchange of a section of 
the micellar network between junctions. Time-scales Trs or r^i can then be defined as the 
lifetime of a typical network strand before destruction by such a process. In the reversible 
scission case (II. lip still holds, now with L the mean strand length in the network [15]. 

In the presence of rings, the three reaction schemes of Figure 11.21 remain applicable in 
principle. It is then notable that the chain number J2nc{N), though not the ring number 
J^N^riN), is still conserved by the two interchange processes. Whenever open chains are 
present, reversible scission is needed for them to reach full thermal equilibrium [23]. 



1.2.6 Parameter Variations 



As stated previously, the static mean-field theories given above (in Sections 11.2.11 - 11.2.31) 
take as their parameters E, Ip, oq, Vq. Also relevant is the excluded volume parameter w [5,6]. 
This controls the strength of repulsions between sections of micelle; for hard core interactions 
this is a function of Ip and Oq, but in general w also depends on all local interaction forces 
between sections of micelles. None theless, with in mean-field, this parameter only affects the 
purely 0-dependent term -Fo(0) in (11.1111.5(11.61) and hence has no effect on the mean micellar 
length L or the size distribution c(L). (The most important role of w is, in fact, to control 
the crossover to the scaling results discussed in Section [1.2.41 above.) All of the parameters 
E, Ip, oq, w in principle can have explicit dependence on the volume fraction 0. This certainly 
occurs in ionic micellar systems at low added salt, where the ionic strength depends strongly 
on itself and modulates directly parameters such as E and Ip. Ion-binding and similar 
effects can also be strongly temperature dependent. Similar remarks apply to the reaction 
rate constants krs, kei considered in Section 11.2.51 above, and hence also to their activation 
energies Ea = —dl nk/djS . The rheological consequences of these parameter variations are 
discussed in Section 11.4.41 
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1.3 Theoretical Rheology 



Since microscopic models for giant micelle rheology draw strongly from earlier progress in 
modelling conventional polymers, we review that progress briefly here. (See [6] for a definitive 
account.) In doing so we can also establish some of the concepts and terms used in rheology 
- a field which remains regrettably foreign to the majority of physics graduates. 



1.3.1 Basic Ideas 



Rheology is the measurement and prediction of flow behaviour. The basic experimental 
tool is a rheometer - a device for applying a controlled stress to a sample and measuring 
its deformation, or vice versa. However, in recent years a variety of rheophysical probes, 
which allow simultaneous microscopic characterisation or imaging, have been developed [34, 
35]. For the complex flows that can arise in giant micelles, these enhanced probes offer 
important additional information about how microstructure and deformation interact. Many 
rheometers use a Couette cell, comprising two concentric cylinders, of radius r and r + h 
with the inner one rotating. (See Figur e 11.191 below for an illustration of this geometry.) 
Others use a cone-plate cell (Figure 11.30 where a rotating cone contacts a stationary plate 
at its apex, with opening angle 6. In the limit of small h/r or small 6, each device results 
in a uniform stress in steady state; in each case, the shear stress can be measured from the 
torque. Some cone-and plate devices can also measure 'normal stress differences' defined 
below. 




Figure 1.3: A cone-plate rheometer. The sample (black) sits between a rotating cone (white) 
and a solid plate (grey). 



Statistical Mechanics of Stress 



We shall use suffix notation, with roman indices and the usual summation convention, for 
vectors and tensors; letters a...w can therefore stand for any of the three cartesian directions 
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X, y, z. Greek indices will be reserved for labels of other kinds. 

Consider a surface element of area dA with unit normal vector n^. Denote by dFi the 
force acting on the interior of the surface element caused by what is outside. If is reversed 
(switching the definitions of interior and exterior), then so is Fi\ this accords with Newton's 
third law. Writing the usual vectorial area element as dSi = riidA, we have 

dFi = CTijdSj (1-13) 

which defines the stress tensor aij. This tensor is symmetric. The hydrostatic pressure 
is defined via the trace of the stress tensor, as p = —era/ 3; what matters in rheology the 
(traceless) 'deviatoric' stress afj^^ = (7^ +p6ij. This includes all shear stresses, and also two 
combinations of the diagonal elements, usually chosen as the two normal stress differences, 

Ni = arcx-<Jyy ; N2 = (Jyy-azz (1-14) 



For simplicity we assume pairwise interactions between particles. (The choice of what we 
define as a particle is clarified later.) The force /"^ exerted by particle a on particle /? then 

depends on their relative coordinate rf^ (measured by convention from a to P). But this pair 
of particles contributes to the force dFi acting across a surface element dSi only if the surface 
divides one particle from the other. The probability of this happening is dSivf^ /V where V 
is the volume of the system. (This is easiest seen for a cubic box of side L with a planar 
dividing surface of area A = L"^ with normal Ui along a symmetry axis; Figure [T3D- The 
separation of the particles normal to the surface is clearly i = rf^Ui, and the probability 
of their lying one either side of it is then just i/L, which can be written as Ar°'^ni/V. 
Accordingly, the total force across a surface element dSi is dFi = ~J2at3dA{rj^nj)f^^/V 
which by definition acts outward (hence the minus sign). Bearing in mind (11.131) . this gives 



cr,: 



-V- 



al3 



:i.i5) 



where the average is taken over pairs and p is the mean particle density. 




Figure 1.4: Contribution from a polymer 'subchain' to the stress tensor. The endpoints of 
the chain can be viewed as two particles, with the chain in between supplying a 'spring force' 
between them. 
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An example is shown in Figure 11.41 where a polymer 'subchain' is shown crossing the 
surface. At a microscopic level, one could choose the individual monomers as the particles, 
and their covalent, van der Waals, and other interactions as the forces in f ll.lSp . (This is 

often done in computer simulation [36,37].) But so long as the force /"^ is suitably redefined 
as an effective, coarse grained force that includes entropic contributions, we can equally well 
consider a polymer chain as a sparse string of 'beads' connected by 'springs'. At this larger 

scale, the interaction force /"^ has a universal and simple dependence on rf^, deriving from 
an 'entropic potential' f/(rj) = {3kBT/2b'^)riri, where 6^ = {riTi). This is a consequence of 
the well-known gaussian distribution law for random walks, of which the polymer, at this 
level of description, is an example. The entropic potential is defined so that the probability 
distribution for the end-to-end vector of the subchain obeys P(rj) oc exp[— [/(rj)/^^^]; this 
form identifies U as the free energy. The force now obeys 

ff = -dU{rf)/drf = -{3kBT/byf (1.16) 

which gives, using fll.151) . the polymeric contribution to the stress tensor: 

^pol _ ^spr 'iksT ^ 



< = -^-^{r.r,) (1.17) 



Here the average is over the probability distribution P(rj) for the end-to-end vectors of our 
polymeric subchains (or springs); Nspr/V is the number of these per unit volume. In polymer 
melts, contributions such as the one we just calculated completely dominate the deviatoric 
stress. In solutions there may also be a significant contribution from local viscous dissipation 
in the solvent. In this case, although a formula such as fll.lSp still holds in principle, it is 
more convenient to work with (11.171) and add a separate solvent contribution directly to the 
stress tensor. For a Newtonian solvent, the additional contribution is af°^ = 7]^°^{Kij + Kji), 
where Kij is the velocity gradient tensor introduced below. 



Strain and Strain Rate 



Consider a uniform, but possibly large, deformation of a material to a strained from an 
unstrained state. The position vector rj of a material point is thereby transformed into r[\ 
the deformation tensor Eij is defined by r- = EijVj. For small deformations, one can write 
this as Eij = 6ij + e^j so that the displacement = — rj obeys Ui = CijVj. Alternatively 
we may write this as eij = VjUi. Consider now a time-dependent strain, for which Vi = iii 
defines the fluid velocity, which depends on the position r^. We define the velocity gradient 
tensor Kij = VjVi = e-ij] this is also sometimes known as the 'rate of strain tensor' or 
'deformation rate tensor'. If we now consider a small strain increment, eij6t, 

r,{t + 6t) = {6ij + eijdt) rj{t) (1.18) 

The left hand side of this is, by definition, Eij{t + 6t)ri{0) where the time-dependent defor- 
mation tensor Eij{t) connects coordinates at time zero with those at time t. Inserting also 
rj{t) = Ejk{t)rk{0) we obtain Eij = KikE^j, or equivalently 

dEij/dt = e,kEk, (1.19) 



An important example is simple shear. Consider a shear rate 7 with flow velocity along 
X and its gradient along y: then Vi = 'jySix- The velocity gradient tensor is Kij = 'jSi^Sjy, 
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that is, Kij is a matrix with 7 in the xy position and all other entries zero. Solving (11.191) 
for arbitrary 7(t) then gives E^j = 6ij +7(t, t')5ixSjy where is defined as the deformation 

tensor connecting vectors at time t to those at time t', and 'y{t,t') = // 'j{t")dt" is the total 
strain between these two times. 



1.3.2 Linear Rheology 



Linear rheology addresses the response of systems to small stresses. Imagine an undeformed 
block of material which is suddenly subjected, at time ti, to a small shear strain 7. Taking the 
displacement along x and its gradient along y, we then have for the resulting deformation 
tensor Eij = Sij + •ySixSjy. Suppose we measure the corresponding stress tensor CTij{t). 
Linearity, combined with time-translational invariance of material properties, requires that 

CTyx = CTxy = G{t - ti)7 (1.20) 

and that all other deviatoric components of cTjj vanish, at linear order in 7, by symmetry. 
(For example, A'^i(7) = A'^i(— 7), which requires Ni = 0(7^).) This defines the linear step- 
strain response function G{t). This function is zero for t < 0; it is discontinuous at t = 0, 
jumping to an initial value which is very large (on a scale set by Gq, defined below). This 
largeness refiects the role of microscopic degrees of freedom; there follows a very rapid decay 
to a more modest level arising from mesoscopic (polymeric) degrees of freedom. In most 
cases this level persists for a while, making it useful to identify it as Gq, the transient elastic 
modulus. (In models that ignore microscopies, one can identify Gq = G{t 0+).) On the 
timescale of mesoscopic relaxations, which are responsible for viscoelasticity, G{t) then falls 
further. 

Now suppose we apply a time-dependent, but small, shear strain •yit). By linearity, we 
can decompose this into a series of infinitesimal steps of magnitude 'j{t')dt'; the response to 
such a step is da^yit) = G(t — t')j(t')dt'. We may sum these incremental responses, giving 

axy{t)= f G{t-t'm')dt' (1.21) 



where, to allow for any displacements that took place before t = 0, we have extended the 
integral into the indefinite past. Hence G{t) is the memory kernel giving the linear stress 
response to an arbitrary shear rate history. This is an example of a constitutive equation. 
However, the constitutive equation for nonlinear fiows is far more complicated. 

In steady shear 7(t) is constant; therefore from (ll.2ip one has cr^jy(t) = 7 /-cxd ~ t')dt' . 
However, the definition of a fiuid's linear viscosity (its 'zero-shear viscosity', 77) is the ratio 
of shear stress to strain rate in a steady measurement when both are small; hence 

/•oo 

7] = / G{t)dt = lim \G*{u;)/iuj] . (1.22) 

JO "^^0 



This is finite so long as G{t) decays to zero faster than 1/t at late times, which is true in all 
viscoelastic liquids (as opposed to solid-like materials), including giant micelles. 



18 



CHAPTER 1. RHEOLOGY OF GIANT MICELLES 



Oscillatory Flow; Linear Creep 



The case of an oscillatory flow is often studie d. We write 7(t) = 7oe*'^* (taking the real part 
whenever appropriate); substituting in fll.2ip gives after trivial manipulation 

a.,(t)=7oe*"*G*(a;) (1.23) 

where G*{uj) = iuj /q°° G{t)e~^'^^dt\ this is called the complex modulus. The complex modulus, 
or 'viscoelastic spectrum', is conventionally written G*{(jj) = G'{uj) + iG"{uj) where the real 
quantities G' and G" are respectively the in-phase or elastic response, and the out-of-phase 
or dissipative response. (These are called the 'storage modulus' and the 'loss modulus' 
respectively.) Many polymeric fluids exhibit a 'longest relaxation time' r in the sense that 
for large enough t, the relaxation modulus G{t) falls off asymptotically like exp[— t/r]. In this 
case one has at low frequencies G' ~ oj"^ and G" ~ oj. For polymer melts and concentrated 
solutions, as frequency is raised G' passes through a plateau whereas G" st arts to fall; 
eventually at high frequencies both rise again. This is sketched in Figure 11.51 where, as is 
common practice, a double logarithmic scale is used to plot the viscoelastic spectra. 



InG' 
InG" 




In (ji 



Figure 1.5: Artist's impression of the viscoelastic spectrum for a typical polymeric material; 
the storage and loss moduli G' {uj) , G" {u) are solid and dotted lines respectively. 

One can also study the steady-state flow response to an oscillatory stress. This defines a 
frequency-dependent complex compliance J*{oj)-^ however, within the linear response regime 
this is just the reciprocal of G*{uj). Suppose, instead of applying a step strain as was used 
to define G{t) in (11.201) . we apply a small step in shear stress of magnitude (Jq and measure 
the strain response 7(t). This defines a compliance function 'j{t)/ao = J{t) which is the 
functional inverse of G{t) (that is, / J(t)G(t — t')dt = S(t')). To see this, one can repeat the 
derivation of (I1.23P with stress and strain interchanged, to find that J*{u!) = iuj J J(t)e~'^^^dt. 
For a viscoelastic liquid 7(t) rises smoothly from zero, and the system eventually asymptotes 
to a steady flow: 7()f: oo) = ao{t/T] + J^^^). The offset J^^\ measured by extrapolating 
the asymptote back to the origin, is called the steady-state compliance. It can be written as 
J^^^ = tG{t)dt/rf and is therefore more sensitive to the late-time part of G{t) than the 
viscosity rj = /q°° G{t)dt. 
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The Linear Maxwell Model 



The simplest imaginable G(t) takes the form G(t) = Goexp(— t/rjv/) for all t > and is 
called the linear Maxwell model, after its inventor James Clerk Maxwell. Go is a transient 
elastic modulus and tm a relaxation time (in this model, it is the only such time) called the 
Maxwell time. The viscosity is rj = GqTm] note that a Newtonian fluid is recovered by taking 
Go oo and tm ^ at fixed rj. In nature, nothing exists that is quite as simple as the 
Maxwell model: but the low-frequ ency linear viscoelasticity of certain giant micellar systems 
is remarkably close to it (Figure II. 6p . The viscoelastic spectrum of the Maxwell model 
is G *{u) = Gqiujtm / + i^TAi) whose real and imaginary parts closely resemble Figure 
11.61 a symmetric maximum in G" on log-log through which G' passes as it rises towards 
a plateau. This is distinct from ordinary polymers, where the peak is lopsided (with slope 
closer to —1/2 on the high u side), with G' not passing through the maximum (Figure [T75|) . 
Understanding the near-Maxwellian behaviour of giant micelles in linear rheology is one of 
the main achievements of the 'reptation-reaction' models outlined in Section 11.41 below. 
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Figure 1.6: Viscoelastic spectrum for a system of entangled micelles: arguably nature's 
closest approach to the linear Maxwell model, for which the peak in G" is perfectly symmetric 
and G' crosses through this peak at the maximum. Figure reprinted with permission from 
Ref. [38]. 



1.3.3 Linear Viscoelasticity of Polymers: Tube Models 



Figure [LT] shows a flexible polymer. The chain conformation is a random walk; its end-to-end 
vector is gaussian distributed. In both polymeric and micellar systems there are corrections 
to gaussian statistics arising from excluded volume effects at length scales smaller than the 
static correlation length ^. These effects are s creene d out at larger distances [5], and their 
effects in micelles anyway limited (see Section ri.2.4p : we ignore them here. 

Dense polymers are somewhat like an entangled mass of spaghetti, lubricated by Brownian 
motion. The presence of other chains strongly impedes the thermal motion of any particular 
chain. Suppose for a moment that the ends of that chain are held fixed. In this case, the 
effect of the obstacles can be represented as a tube (Figure [L71 A). Because it wraps around 
a random walk, the tube is also a random walk; its number of steps and step-length b 
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Figure 1.7: A polymer chain (light line). Surrounding chains present obstacles that the chain 
cannot cross. These can be modelled by a tube (heavy line). The stress relaxation response 
after step strain is controlled by the fraction of the intial tube still occupied by the chain 
at time t. Frames A-D show the state of the tube at four consecutive times, with vacated 
regions shown dotted. Although the emerging chain creates new tube around itself (not 
shown) this part is assumed to be created in the strained environment, and hence carries no 
shear stress. 



1.3. THEORETICAL RHEOLOGY 



21 



(comparable to the tube diameter) must obey the usual relation (i?^) = A^r^^ where R is 
the end-to-end distance of both the tube and the chain. This distance can be measured by 
scattering with selected labelling, as can, in effect, the tube diameter (or step length), by 
looking at fluctuations in chain position on short enough timescales that the chain ends don't 
move much. However, there is no fundamental theory that can predict 6; in what follows it 
is a parameter. It is quite large, so that chains smaller than a few hundred monomers do 
not feel the tube at all. (The largeness of b remains an active topic of research [39].) In 
what follows we will address strongly entangled materials for which Nt ^ 1, ignoring many 
subtle questions that arise when Nt is of order unity. 

Suppose we now take a dense polymer system and perform a sudden step-strain with shear 
strain 7. The chain will instantaneously deform with the applied strain. Since a deformed 
random walk is not maximally random, but biased, this causes a drop in its configurational 
entropy. Quite rapidly, though, degrees of freedom at short scales (within the tube) can 
relax by Brownian motion. Once this has happened, the only remaining bias is at the scale 
of the tube: the residual entropy change AS* of the chain is effectively that of the tube in 
which it resides. A calculation [6] of the entropy of deformed random walks gives a resulting 
free energy change 



where we identify Go as the transient elastic modulus; this comes out as Go = ^ksTn/b 
where n is the number of tube segments per unit volume. Hence the elastic modulus is close 
to, but not exactly, /c^T per tube segment. 

What happens next? The chain continues to move by Brownian motion, as do its neigh- 
bours. Although the individual constraints may come and go to some extent, the primary 
effect is as if the chain remains hemmed in by its tube (Figure fLTI) . Therefore it can diffuse 
only along the axis of the tube (curvilinear diffusion). The curvilinear diffusion constant Dc 
is inversely proportional to chain length L [5] . Curvilinear diffusion allows a chain to escape 
through the ends of the tube. When it does so, the chain encounters new obstacles and, 
in effect, creates new tube around itself. However, we assume that this new tube, which is 
created at random after the original strain was applied, is undeformed. This turns out to 
be a very good approximation, mainly since h is so large: the deformation at the tube scale 
leads to a local, molecular level alignment that is very small indeed. (Such an alignment 
might 'steer' the emerging chain end so that new tube was correlated with the old; this effect 
is measurable [40], but negligible for our purposes.) This causes the stored free energy AF 
to decay away as AF = G(t)7^/2 where 



Here we identify /i(t) as the fraction of the original tube (crea ted at time zero) which is 
still occupied, by any part of the chain, at time t. (In Figure 11.71 this part of the tube 
is shown with the solid line in each time frame; the remaining, vacated, regions are shown 
dotted.) The problem of finding /i(t) can be recast [5] as the problem of finding the survival 
probability up to time t of a particle of diffusivity D^. which lives on a line segment (0,L), 
with absorbing boundary conditions at each end; the particle is placed at random on the line 
segment at time zero. To understand this, choose a random segment of the initial tube and 
paint it red; then go into a frame where the chain is stationary and the tube is moving. The 
red tube segment, which started at a random place, diffuses relative to the chain and is lost 
when it meets a chain end. This tube segment is our particle, and its survival probability 
defines /i(t). 



AF = -TAS = -Go7^ 



(1.24) 



G(t) = GoAi(t) 



(1.25) 
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It is remarkable that the tube concept simphfies our dynamics from a comphcated many- 
chain problem, first into a one-chain (-|-tube) problem, and then into a one-particle problem. 
The result of this calculation, a good revision exercise in eigenfunction analysis [6], is: 

M^)= E -4-^exp\-nH/rR] (1.26) 
ra=odd 

where tr = L'^n^'^/Dc- This parameter is called the 'reptation time' ('reptate' means to 
move like a snake through long grass), and sets the basic timescale for escape from the tube. 
The calculated fi{t) is dominated by the slowest decaying term - hence it is not that far from 
the Maxwell model, though clearly different from it, and resembles the left part of Figure 
11.51 (To understand the upturn at the right hand side of that figure, one needs to include 
intra-tube modes; see [6].) From this form of n{t) follow several results: for example the 
viscosity is rj = J G{t)dt = GoTrTi'^/12 and the steady-state comphance obeys Jf^Go = 6/5. 
Thus the tube model gives quantitative inter-relations between observable quantities, and 
the number of these relations significantly exceeds the number of free parameters in the 

theory — which can be chosen, in effect, as Go and a diffusivity parameter Dc = D^L. 

The model predicts that rj = GqL^ / {12Dc); since Gq is independent of molecular weight, 
7] at fixed varies as L^ for long chains. The experiments lie closer to ?7 ~ L'^'^, at least 
for modest L, but with a prefactor such that the observations lie below the tube model's 
prediction until extremely large L is attained (at which point, in fact, the data bend over 
towards L^). This viscosity deficit at intermediate chain lengths has, in recent years, been 
successfully accounted for by studying more closely the role of intra-tube fluctuation modes 
and their effects on other chains; see [41]. 



1.3.4 Nonlinear Rheology 



Nonlinear rheology addresses the response of a system to finite or large stresses. In the 
absence of a superposition principle, such as the one that holds for linear response, the 
range of independent measurements is much wider. Nonlinear versions exist of the step- 
strain and step-stress response measurements discussed in Section 11.3.21 and of oscillatory 
measurements in which either stress or strain oscillate sinusoidally (though in the nonlinear 
regime, the induced strain or stress will have a more complicated waveform). 

In nonlinear step strain, a deformation Eij = 6ij + 'jSixSjy is suddenly applied at time ti, 
just as in Section 11.3.21 but now 7 need not be small. Analogous to ( ll.20p we define 

a,y = G{t-ti;^)^ (1.27) 

where a factor of 7 ensures that G(t — ti; 0) = G{t — ti) (so the small-strain limit coincides 
with the linear modulus defined previously). A system is called 'factorable' if G{t — ti; 7) = 
G{t — ti)h{'-f), but this is not the general case. Whereas at linear order all other deviatoric 
components of aij vanished by symmetry, in the nonlinear regime one can expect to measure 
finite normal stress differences Ni,N2, as defined in f ll.Mp . In some cases, including many 
systems containing giant micelles, these quantities greatly exceed the shear stress axy [3]. 

Another key experiment in the nonlinear shear regime is to measure the 'flow curve', that 
is, the relationship (7(7) in steady state. For a Newtonian fluid this is a straight line of slope 
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rj; upward curvature is called shear-thickening and downward curvature shear-thinning. Flow 
curves can also exhibit vertical or horizontal discontinuities: these are usually associated with 
an underlying instability to an inhomogeneous flow, to which we return in Section 11.51 



Nonlinear Step Strain for Polymers 



Imagine a dense polymer system to which a finite strain is suddenly applied. We are thinking 
mainly of shear, but can equally consider an arbitrary strain tensor Eij. As previously 
discussed, the random walk comprising the tube, which describes the slow degrees of freedom, 
becomes non-random. If we define the tube as a string of vectors buf (where a labels the 
tube segment) then the initial u°' are random unit vectors. After deformation 

< Eiju] (1.28) 

where it is a simple matter to prove [6] that the average length of the vector has gone up: 
(lEijU'jDa = X > 1) where the average so defined is over the initial, isotropic distribution. 

The length increment is of order 7^ (for the usual reasons of symmetry; strains 7 and —7 
must be equivalent, macroscopically) but for large strains cannot be neglected. 

This increase in the length of the tube is rapidly relaxed by a 'breathing' motion [6] of the 
free ends (one of the intra-tube modes mentioned previously). This rapid retraction kills off 
a fraction 1 — 1/x of the tube segments, so that in effect n n/x- Retraction also relaxes 
the magnitude, but not the direction, of the mean spr ing fo rce in a tube segment back to the 
equilibrium value. The resulting force according to (11.161) is = {3kBT/b)Eiju'^ /\EijU'^\, 
while the corresponding end-to-end vector of the segment is hEijU^. Substituting these 
results in f ll.l7p gives 

Here the final fi{t — ti) is inserted on the grounds that, after retraction is over, the dynamics 
proceeds exactly as discussed previously for escape of a chain from a tube. 

This stress relaxation is of factorable form (now choosing ti = 0): 

aijit) = ?,nkBTQij{Emn)^-{t) (1.30) 

which defines a tensor Qij as a function of the step deformation Emn- Computing Qij involves 
only angular integrations over a sphere, since the a average in fll.29p is over isotropic unit 
vectors [6]. Expanding the result in 7 for simple shear gives Qxy = 47/I5 + 0(7^); this 
confirms the value of the transient modulus Gq quoted after (11.241) above. In finite amplitude 
shear, Qij is sublinear in deformation: this is called 'strain-softening' and the same physics 
is responsible for shear-thinning in polymers under steady flow. 

There are two ways to explain this sublinearity. One is retraction, leading to loss of tube 
segments. The other is 'over alignment': a randomly oriented ensemble of tube segments 
will, if strained too far, all point along the flow direction. Hen ce none will cross a plane 
transverse to the flow as required to give a shear stress (Figure [L^ . But the second argument 
is fallacious unless retraction also occurs (the number of chains crossing the given plane 
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is otherwise conserved) and indeed crosslinked polymer networks, where retraction cannot 
happen because of permanent connections, do not strain-soften. Like many other predictions 
of the tube model, these ones are quantitative to 10 or 15 percent. Note that the factorability 
stems from the separation of timescales between slow (reptation) modes and the faster ones 
(breathing) causing retraction; close experimental examination shows that the factorisation 
fails at short times. 



Constitutive Equation for Polymers 

Alongside shear thinning, polymeric fluids exhibit several exotic phenomena under strong 
flows; these go by the names of rod-climbing, recoil, the tubeless syphon, etc. [42]. Because 
the behaviour of a viscoelastic material cannot be summarised by a few linear or nonlinear 
tests, the goal of serious theoretical rheology is to obtain for each material studied a consti- 
tutive equation: a functional relationship between the stress at time t and the deformation 
applied at all previous times (or vice versa). The tube model, in its simplest form (involving 
a further simplification called the 'independent alignment approximation, or lAA') has the 
following constitutive equation, due (like so much above) to Doi and Edwards [6]: 



where Qij{Emn) is as defined in fll.30p and denotes the deformation tensor connecting 
the shape of the sample at time t to that at an earlier time t'. (Recall this is found by 
solving (11.191) with initial condition E^^ = 6mn, so it is fully determined by the strain rate 
history.) This is the deformation seen by tube segments that were created at time t'; GoQij 
gives the corresponding stress contribution. The factor fi(t — t') (with nit) obeying (11.261) ) 
is the probability that a tube segment, still alive at time t, was created at the earlier time t' . 

We see that, despite its tensorial complexity, the constitutive equation for the tube model 
(within the lAA approximation, at least) has a relatively simple structure in terms of an 
underlyin g 'bir th and death' dynamics of tube segments. The Doi-Edwards constitutive 
equation (11.311) has formed the basis of a series of further advances in which not only lA A but 
several other simplifications of the tube model have been improved upon - see Section 11.4.31 
and the review by McLeish [39] . Often these more careful theories add no further parameters 
to the model; it is remarkable that, in almost every case, agreement with experiment gets 
better rather than worse when such changes are made. This is a very strong indication 
that the basic concept of the tube mode is very nearly correct - something far from obvious 
when (11.311) was first written down in 1978 [43]. Among its 'unforced triumphs' were Jf"* 
independent of molecular weight; jf^Go a constant not far above unity; zero-shear viscosity 
7] ^ L^ (not far from the experimental exponent); and factorability in step strain with 
roughly the right strain dependence [43]. 



1.3.5 Upper Convected Maxwell Model and Oldroyd B Model 

For some macroscopic purpose s (S ection 11.51 below), constitutive models hke ( ]1.3ip . and the 
analogues presented in Section [TTil for giant micelles, are rather too complicated. Most of the 
macroscopic studies start instead from simpler models which (thanks to various adjustable 




(1.31) 
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parameters) can be tuned to mimic the mi cellar problem to some extent. Some of these 
simpified models can in turn be motivated by the so-called dumb-bell picture, which in fact 
predated the tube model by many years. 

A polymer dumb-bell is defined as two beads connected by a gaussian spring. We forget 
now about entanglements, and represent each polymer by a single dumb-bell, whose end-to- 
end vector is Ri. The force in the spring is = —XRi. (Hence A = SkBT/Nmb'^ where Nm 
is the number of monomers in the underlying chain and bm is the bond length; but this does 
not actually matter once we adopt the dumb-bell picture.) In thermal equilibrium, it follows 

that {RiRj)e = ksTdij/X and we can write f ll.l7p as off" = nDX{RiRj), where hd = N^/V 
is the number of dumb-bells per unit volume. The dumb-bell model assumes that the two 
beads undergo independent diffusion subject to (a) the spring force, and (b) the advection 
of the beads by the fluid in which they are suspended. These ingredients can be combined 

to give a relatively simple equation of motion for (T?°\ as follows. 

First, consider diffusion alone. This would give d{RiRj)/dT = 4:kBT6ij/(. This equa- 
tion says that the separation vector evolves through the sum of two independent diffu- 
sion processes, each of diffusivity D = ksT/^, and hence with combined diffusivity 2D] 
( is the friction factor (or inverse mobility) of a bead. Next, add the spring force: this 
gives a diffusive regression towards the equilibrium value of {RiRj)e = ksTdij/X, that is: 
d{RiRj) / dt = [AksT IQ {6ij — X{RiRj)/kBT). Finally, we allow for advection of beads by 

the flow; on its own this would give i?, = KijRj, from which it follows that d{RiRj) / dt\^g^ = 
Kii{RiRj) + {RiRi)Kji. Combining these elements yields 

J^{R^Rj) = Ku{RiRj) + {RA)Kji + {6ij - X{RiR,)/kBT) (1.32) 

which is equivalent to 

= Kaaf + oI^'k,, + r ^ (G,6,, - af) (1.33) 

where r = C,/ {AksTX) is the relaxation time, and Go = N^X/V is the transient modulus, 
of the system. This is a differ ential constitutive equation, which can also be cast into an 
integral form resembling (11.311) : it is called the 'upper convected Maxwell model' [42]. 

The equations above consider only the polymeric contribution to the str ess. T o this can 
be added a standard, Newtonian contribution from the solvent (see Section 11.3.11 above) 

a,^=af + r]'°\K,,+K^i) (1.34) 

Eq jl.341 defines the so-called Oldroyd B fluid. This model is the most natural extension 
to nonlinear flows of the linear Maxwell model of Section 11.3.21 and so its adoption for 
macroscopic flow modelling in micellar systems, which are nearly Maxwellian in the linear 
regime, is highly appealing. However, this is not enough - in particular it cannot describe the 
spectacular shear-thinning behaviour, and related flow instabilities, seen in these systems. 
The simplest model capable of this is called the Johnson- Segalman model, which will be 
presented in Section 11.51 it reduces to Oldroyd B in a certain limit, but has additional 
parameters allowing a much closer approach to micellar rheology. 

The Oldroyd B fluid is also closely related [42] to the Giesekus model which has sometimes 
been advocated as a versatile modelling tool for macroscopic micellar rheology [44]. Caution 
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is needed however: this can easily become pure curve-fitting i f, fo r instance, the model 
assumes homogeneous uniform flow when (as explored in Section ll.Sp the experimental flow 
curve in fact represents an average over what is an unsteady or inhomogeneous situation. 



1.4 Microscopic Constitutive Modelling of Giant Micelles 

In 1987 one of us [8] proposed an extension of the tube model of polymer viscoelasticity that 
allows incorporation of micellar reactions. This led to a predictive constitutive model for 
viscoelastic surfactant solutions. Here we review the model (Section II. 4. 2p . outline its main 
rheological predictions (Section I1.4.3P and briefly overview the extent to which these have 
been experimentally verified. There follow discussio ns of complexities arising from ionicity 
effects and branching in entangled micelles (Section ll.4.5p . 

Although generally successful in the highly entangled region, the microscopic approach 
initiated by Ref. [8] has not proved easily generalizable to the shear-thickening window 

around the viscoelastic onset threshold 0. Here, one has a system which apparently becomes 
entangled only as a result of structural buildup upon shearing; in some cases there is also 
evidence of nematic or other ordering within the resulting shear-induced viscoelastic struc- 
ture. Attempts to model these phenomena are briefly outlined in Section 11.4.61 Then, we 
address in Section [1.4.81 various 'structural memory' effects, in which material properties of 
micellar systems can evolve on time scales much longer than the Maxwell time. 



1.4.1 Slow Reaction Limit 

Consider first a system of (linear, unbranched) giant micelles for which the kinetic timescales 
Trs, Tei, Tbi govemiug revcrslblc scission and interchange reactions are exceedingly long. After 
waiting this long time, the system will achieve equilibrium with size distribution obeying 
(11. 3p . (Below we assume mean- field theory is appropriate unless otherwise stated.) This 
creates an exponentially polydisperse system, but, if the micelles are entangled, as we from 
now on assume, the system is otherwise equivalent to a set of unbreakable polymers. This 
is because the identity of any individual chain is preserved on the time scale of its stress 
relaxation. Hence, for the purposes of calculating the stress relaxatio n func tion fi(t), defined 
previously by writing G(t) = Go^(t) as the step-strain response in (ll.20p . one has a pure 
polymer problem. Recall that Go, the plateau modulus, depends only on the micellar contour 
length per unit volume; if parameters such as Ip, are held constant. Go ~ 0^ in mean-field 
(or 0^'^ in a scaling approach). 

The simplest approach then is to write the overall stress relaxation function as the length- 
weighted average over all the chains present in the system [8, 11]: 

/.(t) = ^ Lc{L)^iL{t)/ E Lc{L) (1.35) 

L L 

Here /iL(t) is the function defined in (11.260 appropriate to the given chain length L, which 

controls the reptation time in that expression via tt^t/? = L^ jDf.- (Recall that D^, the 
curvilinear diffusivity, is L-independent, though it does depend on 0, /p, oq and the solvent 
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viscosity r] which controls the local drag on a chain.) An estimate of (11.351) gives [8] 



This relaxation function has a characteristic relaxation time given by the reptation time for 
a chain of the average length (we abbreviate = tji{L)) but, in contrast to th e result ( 11.26^ 
for monodisperse chains (let alone the linear Maxwell model of Section 11.3.21) it represents 
an extremely nonexponential decay. 

The above crude result for G{t) can doubtless be much improved by applying modern 
'dynamic dilution' concepts which account for the removal of constraints comprising the 
tube around a long chain, on the time scale of reptation of the shorter ones [39]. (We have 
not seen this worked out for the specific case of exponential polydispersity; but see [45] for an 
indication of the state of the art.) A strongly nonexponential relaxation would nonetheless 
remain; indeed such decays are a well-known experimental signature of polydispersity. The 
experimental observation of a near monoexponential relaxation (e.g.. Figure fIT6l) in many 
viscoelastic micellar systems thus proves the presence of a different relaxation mechanism 
from simple reptation [11]. On the other hand, there are a number of experimental systems 
where results similar to (11.361) are observed [3]; this can be taken as evidence that t/j(L) ^ 
Trsi Tei, Tfoj, SO that miccllar reactions have negligible direct effect on stress relaxation, in those 
systems. 



1.4.2 Reptation-Reaction Model 

We now define r?,, the mean breaking time for a micelle, as the lesser of r^s in (11.111) for 
reversible scission, or Tei in (I1.12p for end interchange. (Bond interchange is dealt with 
separately below.) This is the lifetime of a chain before breaking, and also, to within a 
factor 2 in the case of end-interchange, the lifetime of an end before recombination. We have 
assumed that whichever reaction is faster, dominates. (It would be a slight improvement to 
add both channels in parallel but, given that they will have different activation energies, one 
probably has to fine-tune the system to make them have comparable rates.) We also define 
a parameter ( = Th/rpi, the ratio of breaking and reptation times. When this is large, one 
recovers the results of the preceding section. 

Ref. [8] proposes that, when ( is small, the dominant mode of stress relaxation is as 
shown in Figu re 11.81 The stress relaxation function /i(t) is, just as for unbreakable chains 
(Section ll.3.3p . the probability that a randomly-chosen tube segment, present at time zero, 
survives to time t without a chain end passing through it. However, the original chain ends 
do not survive long enough for ordinary reptation to occur; instead, each tube segment has 
to wait for a break to occur close enough to it, that the new chain end can pass through the 
given tube segment before disappearing again. The distance / an end can move by reptation 
during its lifetime ti, obeys Dc{L)P ~ Tb, hence (l/L)"^ ~ Th/rji where = tr{L) as defined 
previously. The waiting time r for a new end to appear within / is t^L/I. This gives, for 
C <^ 1, a characteristic stress relaxation time 



which is the geometric mean of the timescales for breaking and for (unbreakable) reptation. 
Moreover, if we also define C = Tb/T-, then in the limit C ^ 1 of rapid breaking, both 



/^(t)^exp[-A(t/r^)V4] 



(1.36) 



(1.37) 
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Figure 1.8: The stress relaxation process in the reptation-reaction model with reversible 
scission. For end interchange, the new end is instead created by attack of another chain 
end (which connects to the right hand fragment) and destroyed by its own reaction with the 
central part of another chain. Figure adapted from Ref. [8]. 



the length of the particular micelle seen by a tube segment, and the position of the tube 

segment within that micelle, are randomized by the reaction kinetics of order ( times 
during the stress relaxation process itself. This causes a rapid averaging to occur, so that 
all tube segments experience near-identical probabilities for stress relaxation; there is no 
dispersion in relaxation rates and accordingly, in this limit, the resulting relaxation function 
fi{t) = exp[— t/r] is a purely Maxwellian, mono-exponential decay [8]. 

For modest (, deviations from the Maxwellian form are of course expected; these have 
been studied numerically using a modification of the stochastic diffusion process described 
prior to fll.26p . (Some results for ( of order unity are presented in the next section.) Note 
also that as ( falls below where Nt is the number of tube segments on the average 

chain (a quantity that is often of order 10-50 in micelles) there is a crossover to a new regime. 
In this very rapid breaking regime, the dominant motion of a chain end during its lifetime is 
not curvilinear diffusion, but a more complex motion called 'breathing' in which the length 
of the chain in its tube has fl uctua tions. This motion is well understood for unbreakable 
chains [6] and gives instead of fll.37p [8] 



r ^ T^^^rN'J' (1.38) 



The deviations from a pure Maxwellian relaxation spectrum in this regime have also been 
studied [50], within a model that allows partial access to the high-frequency regime, in which 
stress relaxation involves not only breathing but the intra-tube 'Rouse modes' [6]. This 
gives a high-frequency turnup in the viscoelastic spectra G'{u!),G"{u!), as already depicted 
schematically for unbreakable polymers in Figure 11.51 Such a turnup does occur i n m icelles 
[50] although it can lie beyond the experimental window (as is the case in Figure [L6l) . 
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Bond Interchange 

As mentioned in Section ll.2.5[ bond-interchange does not directly create or destroy chain 
ends and so is less effective than reversible scission or end-interchange at causing stress 
relaxation. Nonetheless, enhancement of relaxation does occur, because bond interchange 
will occasionally bring what was a distant chain end very close to a given tube segment, 
allowing relaxation to proceed faster than on a chain undergoing no reactions. This requires 
the interchange event to create a chain end no further away than I (which was defined as the 
curvilinear distance a chain can move during its lifetime); the waiting time for this is not 
Tb = {kbiL(j))~^ but TbL/l. As shown in [32] the result is to replace (I1.37P with 

r = r'/V,/' (1.39) 

There is a second effect, discussed already in Section [1.2.51 which is the 'evaporation' of the 
tube caused by those bond interchange processes whose effect is to pass one chain through 
another. Closer inspection shows that this does not affect the regime just described in which 
chain ends are reptating on the time-scale of their survival, but alters the analogue of (11.381) 
where they move by breathing modes on this time scale. For details, see [32]. 



Constitutive Equation for Giant Micelles 

We now turn from the stress relaxation function fi{t) (which, alongside Go, is enough to 
determine all linear viscoelastic properties) to the nonlinear constitutive equation of the 
reptation-reaction model. This was first worked out in [46]. We assume A^^^ < C ^ 1 so 
that the linear response behaviour is Maxwellian with relaxation time obeying (ll.37p . and, 
more importantly, all tube segments are governed by the same relaxation dynamics. We also 
assume that the rates of micellar reactions are unperturbed by shear; in the highly entan- 
gled regime, interesting rheology arises even at modest shear rates, so th is is a reasonable 
approximation. (We will revisit it for barely-entangled systems in Section [1.4.61 ) 

As calculated in Ref. [46], the constitutive equation for giant micelles is written in terms 
of the deviatoric part of the polymer stress as 

af''^- = ^Go[W^.,-5.,/3] (1.40) 



W,,it) = f Bivit'))exp\- fv{v{t")) 

J —oo L Jt' 



Q,,(El)cit' (1.41) 



vit) = W,,it)K,^it) (1.42) 
Q.AEmn) = (1.43) 



Eirn.Ur. 







In f ll.40p . Wij{t) = {uiUj) is the second moment of the distribution at time t of unit vector 
orientations Ui for tube segments; this controls the polymeric stress in a manner similar 
to that discussed in Section 11.3.41 f or un breakable chains. (Indeed, this equation represents 
essentially the deviatoric part of fll.l7p .) The central constitutive equation is (11.411) . in 
which B and V are birth and death rates for tube segments. These rates obey complicated 
equations [46] but are well approximated for v > by T> = l/r + v, B = 1/r; and for f < 
hj V = l/r, B = 1/t — V. Here v(t), which obeys (ll.42p . is the rate of destruction of tube 
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segme nts by the same retraction process as was outlined for unbreakable chains in Section 
ll.3.4j this is positive in most flows but can become negative if a flow is suddenly reversed [6], 
in which case chains instead spill out of their tubes causing an addition to the birth rate B. 
In the absence of flow (retraction rate f = 0), the birth and death rates are of course equal 
(by detailed balance), and both are given by 1/r, which is the death rate of a tube segment 
by the reptation-reaction process shown in Figure II. 81 

The physical content of (11.411) is that the stress in the system at any time t is found by 
integrating over past times t' the creation rate B{v{t')) for tube segments, multiplied by an 

exponential factor which is the survival probablity of these segments up to time t, times Qij 
which is the stress contribution of such a surviving segment (allowing for both its elongation 
and orientation by the intervening deformation -E^'„) . The quantity Qij is in turn calculated 
in (11.431) where the average is over an isotropic distribution of initial tangent vect ors Uj . 
This is a close relative (but not identical) to Qij defined for unbreakable chains in (ll.3ip : 
the difference lies essentially in the independent alignment approximation, lAA, which we 
do not require here [46]. We see that, in the limit of rapid breaking C ^ 1, the coupling of 
reactions into the reptation mechani sm pr oduces a constitutive equation that is hardly more 
complicated than the standard one (11.311) for entangled polymers and, because lAA is not 
required, achieves this with slightly fewer approximations than in the case of unbreakable 
chains. Moreover there is no worry about polydispersity (which plagues comparison of (11.311) 
with experiment) since this is erased by the fast averaging process. Giant micelles in the 
fast-reaction regime therefore offer a very interesting arena for testing quantitatively the 
fundamental ingredients of the tube model and related concepts in polymer dynamics [46]. 

Note that the reptation-reaction model, as embodied in (ll.40f[L^3|) . describes entangle- 
ments at the level of 'first generation' tube models (as surveyed in [6]). Since these equa- 
tions were first proposed for micelles in 1990, there has been significant work on developing 
'second generation' tube models which give predictions for unbreakable polymers that are 
significantly closer to experimental observation s. Th e relevance of these 'second generation' 
tube models to micelles is discussed in Section 11.4.31 

In t he lin ear viscoelastic limit, Eqs. (ll.40f[L^3l) reduce to the linear Maxwell model, which 
obeys ( 11.2ip with G(t) = Goexp[— t/r]. Accordingly all lin ear viscoelastic quantities reduce 
to those of a pure Maxwellian fluid with r obeying fll.37p. However, the full constitutive 



model is, in both structure and content, quite unlike (I1.33P for the upper convected Maxwell 
model, which would be — to an empirical rheologist — a 'natural' nonlinear extension of the 
linear Maxwell model. (There are other nonlinear extensions which do show such effects, 
including the so-called 'co-rotational' Maxwell model. However the way stress elements move 
with the fluid in these models is not natural for polymers [42].) 



1.4.3 Primary Rheological Predictions 



As discussed previously, we distinguish primary predictions, which do not require much input 
in terms of how parameters such as E, krs, Ip depend on concentration or temperature, from 
secondary predictions, which do require this information. (For the latter, see Section [1.4.41 ) 
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Linear Spectra: Cole-Cole Plots 



The reptation-reaction model predicts that, in the rapid breaking regime, ii{t) approaches 
an exponential form as described above. Such behaviour has by now been reported dozens 
of times in the literature (see e.g. [3,38,47,48], and many more recent papers). Its prediction 
from a microscopic model is nontrivial: while many simple models, such as that of Maxwell 
himself, 'predict' mono-exponential decay, this is actually an assumption, not a result, of 
such models. In contrast, the reptation-reaction mechanism quantitatively explains the 
monoexponential behaviour, in terms of the rapid averaging of tube-segments over chain 
length and position within the chain on timescale = (t. 

By numerical methods [8,49,50], the model can also predict systematic deviations from 
the Maxwellian form for ( of order unity. In an experimental system where the rapid breaking 
regime does not cover the entire phase diagram but can be left by varying a parameter (salt 
concentration, for example), these predictions offer additional tests of the model, particularly 
in cases where can be estimated independently [31]. Such deviations are rendered most 
visible in the so-called Cole-Cole representation [31,49] whereby G"{uj) is plotted against 
G"{u)] for exponential fi{t) the result should be a perfect semicircle. For systems well within 
the rapid-breaking regime, this has been confirmed with remarkable accuracy in a large 
number of cases (see, e.g., [11,31,33,38,48,51,52]). Significant deviations are predicted [49] 
at ( > 0.4. In at least one case, where Tb is independently determined by temperature jump 
on the same system, a fit between the experimental and theoretical Cole-Cole plots gives good 
agreement on the value of C = rb/r [31]. In more recent measurements, the breaking times 

determined by T-jump at volume fractions below and around (j) were found to extrapolate 
with good coherence to those measured by rheology at higher concentrations [53] ; see Figure 
11.91 This suggests, not only that the reptation-reaction model is correct for this class of 
systems, but that breakage in it takes place by reversible scission, not end-interchange (see 
Section [1.2.51 above). In some systems {e.g., [51]) deviations from the semicircle involve a 
sharp upturn at high uj; this can be understood qualitatively in terms of the intervention of 
intra-tube 'Rouse modes' [50]. For systems where the deviation from the semicircle has a 
region of negative slope prior to the upturn (as in Figure [L9|l . the location of the minimum 
can be used [50] to estimate the mean micellar chain length L. This was done for the system 
of [31], and a physically reasonable trend for L{(f)) found [50]. 

In summary, the linear viscoelastic behaviour of entangled micellar systems, across a 
wide range of different chemical types (see, e.g., [47]) shows a regime of strongly Maxwellian 
relaxation in accord with the reptation-reaction model. The leading shape corrections to the 
relaxation spectrum as one departs this regime are also well accounted for, in many systems, 
by that model. These statements involve no knowledge about how the kinetic and structural 
inputs to the model vary with concentration, ionic strength, or temperature. 



Nonlinear Rheology: Shear-Banding 



The micellar constitutive equation f ll.4ip can be solved for nonlinear step strain; results show 
generally similar features to unbreakable polymer rheology [46]. In practice comparison to 
these predictions are often complicated by strain-induced instabilities an d/or s train-history 
effects [3] , the latter perhaps related to those discussed below in Section 11.4.81 The overall 
similarity between the micellar constitutive equation and that for unbreakable polymers 
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Figure 1.9: Upper panel: comparison of Trs measured by T-jump and rheology in a 
CePyCl/NaSal (2:1) system on a log plot. The horizontal axis is log concentration with 
data spanning 0.001-0.07 g/cm^. Lower panel: fit of Cole-Cole data, viz G"{uj) vs. G'{(jj) on 
a parametric plot (open symbols) to the numerical model of [49] (dots). Data is normalized 
by the radius of the osculating semicircle. Note that the turnup, caused by Rouse modes [50] 
is not part of the fitted data. Figures adapted from [54], courtesy C. Oelschlaeger. 
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Figure 1.10: Flow curves for reptation-reaction model: solid line, by solution of fll.40l - [L^3|) : 
dashed line, with additional quasi-newtonian stress calculated (with one fit parameter) as 
per [66]; dotted-line, top-jumping shear-banded solution (unaffected by that parameter). 
Figure courtesy N. Spenley [61]. 

nonetheless means that most of the nonlinear viscoelastic functions (normal stress coefficients 
\l/i,2 = lini7->o(A'^i,2/7^), etc.), should be broadly similar to those for unbreakable chains [46]. 

The main arena for comparing nonlinear predictions of the microscopic model with ex- 
periments on micelle s ha s involved 'steady' flow. This flow is not always literally steady, 
however; see Section II. 5[ Steady flow, at least if homogeneous, is fully characterized by the 
'flow curve' (7(7), which relates shear stress to strain rate in steady state, and the normal 
stress difference curves A'^i(7) and A^2(7)- Until the early 1990s, studies of the flow curve for 
giant micelles simply assumed homogeneity of the flow. Checking for this has since become 
much easier with a variety of modern techniques [34,35,55-58]. 

A very striking observation, reported first in a CPySal/NaSal system (compare Figure 

II. 6p [3,38], was that above a certain strain rate 7^, the shear stress a attains a plateau 
value 0" = dp, remaining at this level for at least two decades in 7 > 7p. At the same 
time the normal stress difference A^^i continues to increase almost linearly. This represents 
shear thinning of a quite drastic kind. It is exploited in technologies such as hand lotions 
and shampoos, allowing a highly viscous liquid to be pumped or squeezed out of the bottle 
through a narrow nozzle; in a non-shear-thinning fluid of equal viscosity the bottle would be 
likely to break first [59]. Such behaviour in strong shear flows is as ubiquitous a feature of 
the nonlinear rheology of giant micelles as is the Maxwellian spectrum in the line ar rheology . 
Its explanation came in [60] , where t he rep tation-reaction constitutive equations fll.40f[L^3!) 
were solved in steady state (Figure [1.101) . It was found that the shear stress a has, as a 
function of 7, a maximum at (7, cr) = (2.6/r, 0.67G'o)- Such a nonmonotonic flow curve is 
known to be unstable [60-62]; but a steady flow can be recovered by developing 'shear-bands'. 

For a shear-thinning system such as this, the shear bands comprise layers of fluid with 
unequal strain rates but equal stress, their layer normals in the velocity gradient direction 
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Figure 1.11: Comparison of reptation-reaction prediction for shear banded flows [60] with 
observations of [3]. Dots are experiments for shear stress a (below at right) and normal 
stress difference Ni. Solid curve is shear stress in the reptation-reaction model assuming 
top-jumping but with no adjust able p arameters. The dashed curve for normal stress has one 
fitting parameter (as per Figure [TTTUI) . Figure courtesy N. Spenley [61]. 



[60]. In this way, the decreasing part of the flow curve is bypassed by coexistence between 
two bands, one at low 7 = 71 and one at high 72, each of which is on an increasing (hence 
normally stable) part of the curve. The relative amounts of the two bands arrange themselves 
to match the macroscopically averaged imposed strain rate. Assuming the nature of the 
coexisting states does not vary as th eir a mounts change (an assumption that ignores coupling 
to concentration fields, see Section [l.Sp . this gives a horizontal stress plateau as observed 
in [3]. The value of the plateau in CpySal/NaSal could be reproduced by assuming 'top 
jumping', in which one of the coexisting states is at the maximum of the flow curve described 
previously (giving ap = 0.676*0) • However, it is now known that the mechanism for selecting 
coexisting stresses is more complicated than this [63]; top-jumping gives an upper estimate 
for (Tp, and lower values have since been seen in many other materials {e.g., [64,65]). 

The theory of [60] assumes that the high shear branch (which in the system of [3] is at 
72 > lOOOs"^) is not purely Newtonian, as it would be if the micelles were smashed to pieces. 
Instead this phase is taken to contain highly aligned micelles, which can contribute a very 
large Ni while maintaining a small shear stress a. (This requires an improve ment to the 
tube model at high stresses [66], giving a linear (7(7) at high shear rates. Figure [T.IOI ) This 
theory accounts for the observation in [3] that Ni continues to increase almost linearly with 
7 throughout the shear-banding plateau at constant (7(7) = ap (Figure [1.111) . 

Since its first prediction in [60], the evidence for shear-banding in viscoelastic micellar so- 
lutions has become overwhelming [34,55-58,67]. Some authors have preferred interpretations 
involving a shear-induced phase transition to a nematic state [65,67]. However, this is not 
excluded by the above arguments, which make negotiable assumptions about the physics of 
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the high-shear branch at 7 = 72. Such 'semi-thermodynamic' mechanisms of shear-banding 
can also explain of the 'unfolding' of the nonmonotonic flow curve on varying or other con- 
ditions, whereby the shear-banding region gets narrower and eventually disappears [65]; but 
a crossover from the Maxwellian rapid-breaking regime into the unbreakable polymer limit 
can explain this too (see below). Even more complicated identifications of the high shear 
branch are not ruled out; these might include a lon g-lived gel phase similar to that induced 
by flow in the shear-thickening region (see Section ri.4.6p . In all cases, there can also be a 
concentration differential between the two bands. This is not addressed in the rept at ion- 
react ion model, and causes a ramp [68] rather than a plateau, in shear stress (Sections 11.5.21 
ll.5.3p . Such couplings to nematic and/or concentration fields take us beyond the limit of 
practical calculability within such a model; these phenomena can instead be addressed with 
the macroscopic approaches of Section 11.51 

Note again that the 'natural' nonlin ear ge neralization of the linear Maxwell model, namely 
the upper convected Maxwell model of fll.33p . does not predict shear banding. The widespread 
observation of this phenomenon [3,55-58] is thus not implicit in the linear response of the 
system, but confirms a quite separate primary prediction of the rept ation- reaction model. 



Second Generation Tube Models: Convective Constraint Release 

One might ask whether such banding phenomena can be expected in ordinary, unbreakable 
polymers. Though seemingly robust for micelles, there is only very limited (and contro- 
versial) evidence for shear banding in polymers [69]; the consensus view is that, even for 
monodisperse chains, shear banding does not occur in the unbreakable case. Certainly there 
is no experimental evidence for shear banding in even moderately polydisperse polymers (let 
alone the exponential polydispersity represented by (11. 3p ). Since this physics should be 
recovered in the unbreakable limit, the reptation-reaction model includes among its robust 
primary predictions that the shear banding instability, seen in the Maxwellian regime, should 
disappear as one crosses into the highly non-Maxwellian regime a.t ( = r^/r/j ^ 1. 

Intriguingly however, the 'fi rst generat ion' tube models, which underly the reptation- 
reaction constitutive equations (ll.40f[L43ll . do actually predict a nonmonotonic flow curve 
for monodisperse unbreakable chains in steady shear. This has formed a large part of the 
inspiration for the development of 'second generation' tube-model constitutive equations. 
These deploy the concept of 'convective constraint release' (CCR) whereby the nonlinear 
stretching and retraction of chains in strong flows reduces the strength of entanglements 
[39,70,71]. Most of these models entail at least one phenomenological parameter describing 
the strength of the CCR effect; in [71], for example, the relevant parameter is called c,^. 

It is important to ask whether the CCR concept, if consistently applied to the case of 
entangled micelles, restores monotonicity of the flow curve in that case also. If it did, this 
would create a direct conflict between the modelling of unbreakable chains (which require 
CCR to restore monotonicity and avoid shear bands) and micelles (which require nonmono- 
tonicity to explain shear bands). This issue is addressed directly in [71]. It is found that, 
on varying the parameter c^, (which vanishes in the first-generation theories) a significant 
range exists where the flow curve for unbreakable chains is monotonic but that for micelles 
is not. The requirement that the model admit shear-banding in micelles is then used by the 
authors of [71] to limit c^, to that window. 

The second-generation constitutive model for micelles derived in [71] directly combines 
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Table 1.1: The exponent g (where r ~ 0^ and rj ~ 0^+^-^) for various regimes. 



Reaction 



C > 1 1 > C > 



reversible scission 
end interchange 
bond interchange 



3^ 
3.4 
3.4 



1.4 
1.2 
1.7 



0.9 
0.6 
0.3 



the physics of CCR with the reptation-reaction picture that led to (11.4011021) . The resulting 
equation s et is less transparent, however; we refer the reader to [71] for details. In brief, 
retaining (11.401) . these authors derive a differential equation for a quantity fij{x) whose 
integral over x plays the role of Wij in (ll.40f[L42|) . This differential equation involves a 
convection term related to (ll.42p . a stretching/orientation term related to (11.431) . and a 
new CCR term (oc c^) which has a relatively elaborate dependence on the flow history. 
Presu mably, although it is not proved explicitly in [71], this formulation reduces to (11.401 - 
11.421) when the CCR term is switched off. 

Although a careful numerical comparison has not yet been made for anything other than 
steady shear, it seems likely that for modest values of Cjy (chosen to retain the nonmono- 
tonicity of the flow curve) the primary rheological predictions of this second-generation 
micel lar model b roadly coincide with those found from the original reptation-reaction model 
Eqs. (ll.40f[L^3l) . This expectation is based on more extensive comparisons that have been 
done between first- and second-generation models in the case of unbreakable chains [71], 
where rheological features unconnected with flow-curve monotonicity are not strongly al- 
tered. 



1.4.4 Secondary Predictions: Concentration and Temperature Dependence 

If E, krs (or kei, kbi, as applicable) and Ip are independent of concentration, a number of 
secondary rheological predictions can be obtained for the concentration dependence of r (and 
hence viscosity r] = Gqt) and related quantities. However, these depend on: (a) whether 
the scaling or mean-field theory applies; (b) whether reversible-scission, end-interchange, 
or bond- inter chang e dominates; (c) whether one is in the pure reptation-reaction regime 
governed by (11.371) . or in the breathing-reaction regime obeying (11.381) . Accordingly such 
laws for concentration trends are at best a rough guide, from which only rather strong 
deviations can be taken as firm evidence against the underlying model. 

Nonethless, we present in Table 11.4.41 'best estimate' scaling laws for an exponent g 
relating the Maxwell time to concentration (r ~ 0^) for different reaction schemes in the 
three different regimes governed by unbreakable chains, reptation-reaction, and breathing- 
reaction processes. These best estim ates, r eproduced from [32] , make use of scaling exponents 
in the absence of rings (see Section 11.2.41 above) for both static properties and the effective 
tube diameter. (In the scaling perspective, the latter is the mesh size, ^ [5].) However, the 
results would not be very different if mean-field estimates were applied consistently instead. 
To convert from r to the viscosity, one multiplies by Go ~ 0^'^- 



We note that g is always positive; in no regime does r fall with 0. This is reasonably 
robust, but could be altered if for some reason the reaction rate for the relevant reaction 
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scheme {krs, say) was strongly increasing with 0. Observation of small positive exponents is 
not disproof of the model, but might be evidence that a system is in the breathing regime 
with interchange kinetics. Thirdly, even if kinetic parameters such as krs,ei,u and static ones 
such as -E, /p are 0-independent as assumed, and breaking is rapid, there can be crossovers 
between six different regimes (Table [T.4.41) : these might normally prevent the experimental 
verification of simple power laws of any kind. Even more complex trends ca n be ra tionalized, 
a posteriori, once ionicity effects are expressly considered [47] ; see Section 11.4.51 Although 
perhaps more informative in being broken than obeyed, the scaling laws in Table 11.4.41 do 
give reasonable agreement with some of the simpler experimental systems (e.g., [11,51,72]). 

One can make further interesting predictions by assuming that E is independent of T 
and also that the activation energy E^ (for whichever reaction is dominant) is likewise T- 
independent. Note that Ea and E are not the same, even for reversible scission (unless there 
happens to be no activation barrier for fusion of chain ends). Nonetheless, since E and Ea 
are likely to be much bigger than the Arrhenius energy Es for the solvent viscosity (which 

controls the curvilinear friction parameter Dc), Equation (11.371) for reversible scission (say) 
predicts that r oc exp[/?£"] with E' ~ {Ea + E)/2. Since Ea is measurable in temperature 
jump, this allows a check on the equilibrium scission energy E:, reasonable values (around 
20k bT) have been obtained in this way, as have reasonable trends for dependence on ionic 
strength and other factors [11,12,53]. In one recent study [12], strong evidence is given for 
the switching off of the reversible scission process at low ionic strength, at least in the specific 
classes of ionic surfactants studied there. This can be rationalized in terms of the electrostatic 
and bending contributions to the activation barriers for micellar reactions. (See [73] for a 
recent use of such ideas in a different context.) At low salt, the micellar breaking time Tf, 
is then controlled by an interchange process, invisible to T-jump, and the activation energy 
comparison with T-jump data described above is inapplicable. 



1.4.5 Role of Branching: Ionicity Effects 

In ionic micellar systems, and also zwitterionic ones, the overall ionic strength and/or the 
degree of specific counterion binding can strongly influence E, Ea, Ip and other parameters 
that were assum ed con stant in the above analysis. While naive apphcation of scalings such 
as those in Table 11.4.41 is thereby precluded, it is possible in many cases to make sense of the 
trends, at least a posteriori, within the physical precepts of the reptation-reaction picture, 
so long as branching is allowed for. Much work along these lines is reviewed in [47]. 

As described for equilibrium statistics in Section [1.2.21 the presence of branching in micel- 
lar chains can alter both static and rheological predictions for their properties. One crucial 
observation is that micellar branch-points are labile; they are always free to slide along the 
length of a micelle. This is quite different from crosslinks in conventional polymers and, 
counterintuitively, means that branching generically decreases the viscosity of the system, 
rather than increasing it. This was first recognized in an important paper by Lequeux [15], 
who showed that, if branch points are present in an entangled micellar network, the curvilin- 
ear diffusion constant of a chain end is = L)c/L, where L now denotes (roughly speaking) 
the distance to the nearest other chain end or junction point in the network. This is be- 
cause a branched network, at scales above the distance between junctions, offers an efficient 
reservoir for micellar contour length. (This can be stored or borrowed by sharing it with 
other network strands, rather than having to move it all the way to the other end of a given, 
unbranched chain.) Accordingly some of the main rheological predictions of the reptation- 
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reaction model can be retained, so long as L carries this new interpretation. Once branching 
is widespread, as shown in Section 11.2.21 one then has L ~ (within a mean-field pic- 

ture); for rheological purposes, the system behaves as z/ micelles were becoming shorter with 
concentration. This can lead to a rheological relaxation time exponent g that, for reversible 
scission reactions, is barely larger than zero (r ~ 0^ = (p^'^^ in a scaling picture [53]). For 
interchange reactions g can even be somewhat negative, in contrast with the results of Table 

MM 

In a careful recent study [53,54], Oelschlaeger et al correlate surfactant hydrophobicity, 
salt, and counterion binding efficiency with the observed g values. These authors find ev- 
idence that the unbranched reptation-reaction model is applicable at relatively low ionic 
strength, high hydrophobicity, and low counterion binding efficiency. They then argue that 
the model of Lequeux [15] explains the falling g values seen when ionic strength is increased, 
surfactant hydrophobicity decreased, or counterion binding increased [47]. Each of these can 
be argued to favour the branching of micelles either through curvature or electrostatic energy 
effects [53]. Moving such parameters further in the same direction can lead to significantly 
negative g and/or a complete coll apse of viscoelasticity [17]. The latter may signify the onset 
of network saturation (see Section ri.2.2p : by this point all potential entanglements have been 
replaced by labile branch points, and resistance to shear is very low. 



1.4.6 Shear Thickening 



Of equal fascination to the Maxwellian and shear-thinning behaviour described above, is 
the phenomenon of shear thickening [12,74-78]. This is seen in a window of volume fraction 

around the onset of viscoelasticity at ~ 0; an initially inviscid or barely viscoelastic system 
is found, after a period of prolonged shearing above a critical shear rate %, to convert into a 
much more viscous state. In some cases the viscosity increase is modest (a factor between two 
and ten) and the state relaxes quite rapidly to the previous one when shearing ceases [74]. 
In other cases, the new phase is a long-lived gel; shear banding is often implicated in the 
formation of such a gel [77]. The thickened phase contains micelles that are nearly fully 
aligned (the optical extinction angle is close to zero); its formation involves a nontrivial 
latency time which in some cases is geometry-dependent [12]. 

At the time of writing there is no microscopic model capable of explaining these phenom- 
ena at anything like the level of a constitutive equation, nor indeed any consensus among 
theorists as to the mechanisms involved. Model-building attempts have been made from time 
to time, based for example on shear-induced aggregation or polymerization of rodlike mi- 
celles [79]. There are also models that couple a gelation transition to shear bands [80,81], and 
some that ascribe the shear thickening directly to ionic or electrokinetic phenomena [82] . All 
these models have drawbacks; for example, aggregation models normally require very high 
shear rates for onset whereas experimental values are low; thickening is sometimes seen in 
nonionic micelles [75]; etc.. In the next Section we briefiy review a radical approach in which 
shear-thickening is speculatively attributed to the presence of micellar rings [23]. 
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1.4.7 Role of Rings 



As discussed in Section 11.2.31 the expected influence of micellar rings is maximal around 
0™""^, the (sadly unknown!) volume fraction at which, in a system of high enough scission 
energy E, a cascade of rings crosses over to a semidilute solution of open chains. The 
reptation-reaction model assumes strong entanglement and hence requires (p ^ 0™°^; in this 
regime any remaining rings are smaller than the mesh size ^ [19] and hence of rather little 
rheological consequence. Thus the role of rings is limited to the concentrati on ran ge below 
the strongly entangled regime, where a cascade of rings is predicted (Section ll.2.3p . 

In the putative cascade-of-rings phase, governed in mean-fleld by (11.71) and present in 
theory at large enough E, it is an open question [23] whether or not the rings interlink so as 
to form a percolating linked network. For unbreakable ring polymers, formation of such a 
structure would create a permanently elastic solid known as an 'olympian gel' [5]. Although 
in a micellar system the rings would have a flnite breakage time, one still expects drastic 
changes to the flow behaviour under any conditions where signiflcant concatenation can arise, 
as seems possible in the 'cascade of rings' regime. Suppose for now that the rings are indeed 
linked for (p ~ (f)'^"'^. If so, then were micellar kinetics suddenly to be switched off, creating 
an olympian gel, the system would acquire some flnite residual modulus Or (assumed small), 
indeflnitely resisting attempts to impose a steady shear flow. Restoring a flnite delinking 
time Ti, this gel becomes a viscoelastic fluid of viscosity rjr = Grn. It is possible that, if Gr 
is small enough, rir remains comparable to that of the solvent; the sample is only marginally 

viscoelastic and in that sense would be considered to have ~ (and identifled as ~ 0* 
if an assumption of linear chains were made). Notice now that the viscoelastic linked-ring 
fluid does not in fact require complete percolation of linked rings; it only requires that linked 
structures extend far enough that their conflgurational relaxation times exceed r;. (If these 
large structures are few and far between, the better is the assumption of small Gr.) 

If such a linked-ring fluid exists around ~ 0, any shear rate 7 > t^~^ will cause strongly 
nonlinear effects. The orientating effect of elastic strains (of order 7r;) on the linked rings 
will alter the reaction rates for bond-interchange and (where present) other reactions. For 
example, pulling two interlocked rings in opposite directions could promote bond interchange 
at their contact point, increasing the creation rate for a larger ring and shifting the mean 
ring size upwards. This polymerization tendency could well cause shear thickening, as could 
tension-induced chain scission, by pushing the equilibrium in (11.91) towards polymerization. 

This scenario, though speculative, does explain several observed features of the shear- 
thickening process. One is a shear-rate threshold for the transition that is lower than any 
reasonable estimate of the reorientation time for marginally overlapped linear micelles [23]. 
Second, since the observed gel phase is fully aligned, a large total strain must be applied 

1 /2 

before it is fully formed. This implies a latency time at least of order TiNji , and potentially 
much longer if the shear-thickening transition is caused by only a slight rate inequality be- 
tween formation and destruction of larger-than-average micelles [23] . The latter mechanism 
could also give a geometry-dependent latency time as reported experimentally [12]. Third, 
any ring-dominated regime requires suppression of open chains in the quiescent state and 
hence large E. Factors favouring large E include raising counterion lipophilicity and using 
gemini (twin-tailed) surfactants; such factors do seem broadly to cause a reduction in 7c 
and enhancement of the viscosity jump. Finally, this picture can tentatively explain some 
strange memory effects seen around the shear thickening transition, discussed next. 
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1.4.8 Structural Memory Effects 



Structural memory is the presence in a system of internal degrees of freedom, other than 
the stress, which relax on a time scale that is at least comparable to the stress relaxation 
time r itself (and in some cases vastly longer). For example, one can find among micellar 
systems some instances where r itself is at most a few seconds; but the value one measures 
for T in repeat experments depends in a complicated manner on the process history of 
the sample over, say, the preceding 24 hours. It is possible, though not yet certain, that 
structural memory plays a major role in the exotic rheological phenomena, such as rheochaos, 
considered in Section 11.51 below. Though hinted at anecdotally from the earliest days of the 
subject [48,74], it is only very recently that structural memory effects in micellar systems 
have been studied in systematic detail [12,78]. 

Among effects observed in [12, 78] are the following. In some shear-thickening micellar 
systems subjected to steady shearing, there is an initial latency time Tiat (seconds or minutes) 
for the thickening to occur. After this, however, the stress level continues to adjust slowly 
over time scales of order hours or days before finally achieving a steady state. If shearing 
is stopped, the stress relaxes quite rapidly but the memory of having been sheared persists for 
times of order r^: if shearing is resumed within this period, a quite different Tiat is measured. 
Moreover the stress level immediately after latency is closer to the ultimate steady-state 
value, and almost identical to it if the switch off period has been short compared to r^. 
Finally, the latency time can also be raised or lowered by a prior incubation at elevated or 
reduced temperature [12,78], even in a sample that has never been sheared. 

These phenomena point to a robust structural property, perturbed by shear but also by 
temperature, as the carrier of structural memory in micellar systems. One such property 
immediate ly spr ings to mind, namely the micelle size distribution. However, as discussed 
in Section ll.2.5[ this relaxes rapidly to its equilibrium form (in mean field, this is the usual 
c(L) oc exp[— L/L]) when any kind of micellar reaction is present. On the other hand, as 
also discussed in Section 11.2.51 the relaxation of the mean micelle length L is contingent on 
the presence of reversible scission reactions. This is because interchange reactions conserve 
the chain number J2l c(-Z^); and, given the fixed shape o f the d istribution, L can only change 
if the chain number does so. As emphasized in Section 11.2.51 chain number is conserved by 
interchange reactions, even when rings are present; but ring number is not itself conserved. 

The structural memory effects reported in [12,78] close to the shear thickening transition 
can be explained in outline if one assumes that (a) both rings and open chains are present 
and (b) reversible scission reactions are very slow. The slow relaxation time is Tg = Trs, the 
time scale for the chain number to reach equilibrium. Stress relaxation is not slow, since 
the faster rate for end interchange or bond interchange reactions will dominate, allowing 
fast relaxation of all quantities other than L, including stress. A system which is sheared or 
thermally treated for a time long compared to Tj-s will acquire a steady-state chain number 
appropriate to those conditions; if conditions are now changed, it will take a time of order 

to relax to the new value. Among things that can vary with chain number are, of course, 
the latency time in the thickening transition; and also the Maxwell time for stress relaxation 
in the quiescent state. 

To understand how the structural memory time can become so long [12,78], note that the 
reaction rates for the three reaction schemes (reversible scission, end-interchange, and bond- 
interchange) all involve different activation energies. For instance, end-interchange involves 
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passing through a state which contains a threefold junction in place of a single end cap; this 
gives an estimate Ej, — E/2 for the relevant activation energy, where Eb is the energy of the 
branch point. The activation energy for scission, however, cannot be less than E, the energy 
to create two end-caps. (An activation energy cannot be less than the energy rise between 
initial and final states.) Using these estimates, we find Trs/Tei ~ exp[{3E/2 — Eh)/kBT]. 
Hence it is indeed quite possible to have a structural relaxation time Ts = Tj-s (governing 
relaxation of the micellar size distribution) that exceeds by orders of magnitude the dominant 
kinetic timescale involved in stress relaxation. The required chemical conditions are those 
that favour large E (small (p) and modest Et. (If the latter is too small, there will be 

significant branching of micelles at concentrations around 0; this will further complicate 
the rheology, but not affect the basic time scale separation under discussion here.) These 
considerations are broadly consistent with the observed trends in Tg on varying the choice of 
surfactant, ionic strength, and other factors [12,23,78]. 

Parts of above scenario could hold equally well for giant micelles in the fully entangled 
regime, where rings are not important. For a system showing slow reversible scission in T- 
jump [12,53], it would be interesting to look for the effects of nonlinear flow, and also thermal 
pre-treatment, on the Maxwell time. If the nonlinear flow creates a shear banding region, 
this might reveal whether the average chain length in the high-shear band is significantly 
perturbed by the flow. A more radical speculation is that at least in some cases, when a 

shear-induced gel phase that forms below persists after stress is removed (and does not 
retain nematic order [67]), this state is in disequilibrium solely through having an enhanced 

L ( of lifetime Tg = r^s). If so, the shear-induced gel state is just another instance of the 
entangled regime of giant micelles, to which the reptation-reaction model can be applied. 
An open issue would be how L acquires its nonequilibrium value during the induction period 
of the shear-thickening transition. 



1.5 Macroscopic Constitutive Modelling 

So far we have discussed microscopic constitutive modelling, which aims to predict rheology 
from an understanding of the micros copic dynamics of the polymer-like micelles themselves. 
We expressed the stres s tens or (11.171) as a sum over micellar chain segments, and derived the 
constitutive equation fll.4ip for its dynamics. Pursuing this approach further, particularly 
to address nonstationary shear-banded flows, becomes prohibitively complicated. Indeed, 
current microscopic models for micelles lack some of the important physics of these flows, 
such as coupling between flow and concentration or collective orientational fluctuations. To 
make progress, we now turn to phenomenological models which, thanks to various adjustable 
parameters, can be tuned to mimic the micellar problem to some extent. 

Before discussing individual models, we sketch the basic features that are common to all 
of them. The stress cxjj is taken to comprise additive contributions from the polymer-like 
micelles and from a Newtonian solvent (as per (11.341) ): 

ai, = af + r^'-\K,, + K^,) (1.44) 

The relevant instabilities are viscoelastic and not inertial in origin, so we work throughout 
at zero Reynolds number. The stress then obeys the force balance equation 

V,a,, - Vjp = V,af;^ + r]'°^V\ - Vjp = (1.45) 
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Here p is an isotropic pressure, which maintains fluid incompressibihty: 

V^Vi = (1.46) 

The viscoelastic stress af°^ is then written as a function of some underlying microstructural 
quantities, whose identities vary according to the system and regime of interest. Common 
choices include the concentration and molecular deformation Wij of the polymeric compo- 
nent; the orientation tensor Qij in nematics; the micellar length distribution P{L), etc.: 

af = afiW,„^,---) (1.47) 

Among these microstructural variables it is important to distinguish "fast" from "slow" 
variables. The former quickly relax to local steady-state values determined by the latter, 
whereas each "slow" variable requires its own dynamical equation of motion. Formally, the 
slowest variables are the "hydrodynamic" ones, which relax at a vanishing rate cx /c" 
{a > 0) for small wavenumbers — >■ 0, either because a long-wavelength distortion of the 
variable costs very little free energy (broken symmetry), or because the quantity is conserved 
and therefore obeys diffusive dynamics (a = 2). In the formal limit k ^ 0, all other variables 
are fast in comparison to these hydrodynamic modes. In viscoelastic solutions however, 
variables that are not strictly hydrodynamic nonetheless relax very slowly. An example is 
the molecular deformation Wij governed by the Maxwell time r (often seconds or minutes). 
It is therefore essential to add to the list of "slow" variables Wij and any other slow but 
formally non-hydrodynamic quantities [83]. 

As discussed in Sec. 11.4.81 above, many micellar systems show a pronounced structural 
memory, with degrees of freedom that relax on a time scale greater than the intrinsic defor- 
mational relaxation time r. Indeed, there are instances in which r is itself at most a few 
seconds, but with a value that varies between repeat experiments in a complicated (though 
reproducible) manner that depends on the process history over the proceeding day or so. This 
points to an underlying structural property that evolves slowly, on a timescale of hours/days. 
One such candidate is the mean micellar length, with a timescale t^s ^ t; another is the 
micellar concentration. 



Consider, then, a scenario in which the relevant dynamical variables are the micellar 
deformation Wij, the micellar concentration 0, and the mean micellar length L, with respec- 
tive relative relaxation timescales r <^ A^/P <C t^s, for a micellar diffusion coefficient V and 
sample size A. The concentration evolves with conserved dynamics of the form 

= VMVfiiVnVm, Wnm, 4>) (1-48) 

with a mobility M, proportional to V. The derivative Dt=dt + viVi denotes rate of cha nge 
in a fluid element convected with the flow field Vi. The direct counterpart of (11.481) for 
nonconserved quantities such as Wij and L is an equation of the form 

DtWij = -Gij{VnV„,, Wnrn, 4>) (1-49) 

r 

Sometimes, however, this differential structure is replaced by an integral form 

L = f dt'M{t-t')g{t') (1.50) 

where g depends on time t' via all relevant quantities Wij{t'), etc.. Only in the case of an 
exponential memory kernel, M.{t) = exp{—t/T)/T, can (11.501) be recast into the differential 
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form of fll.49p . The integral form of (11.501) emphasises, for long-hved variables such as L, 
the structural memory at time t to all earlier times t' < t. 

Note that (11.471) above assumes the polymeric stress to be instantaneously prescribed 
by any given configuration of the underlying microscopic variables, some of which them- 
selves obey "slow" dynamical equations. This approach retains the basic principle that even 
macroscopic models should u ltimately be motivated by microscopies, however oversimpli- 
fied the evolution equations (ll.48f[T750|) become. A simpler, and very common, approach 
consists of directly prescribing the dynamics of the polymeric stress itself by means of an 
autonomous viscoelastic constitutive equation. Fundamentally, though, this should still be 
thought o f in ter ms of the evolution of an underlying microstructural quantity: for example, 
in Section [1.3.51 the equation (I1.33P for stress evolution is actually an instance of (11.491) . 

The basic structure outlined above encompass es the micr oscop ic mode ls discu ssed earlier. 
For example, the reptation-reaction model has (ll.40p for (I1.47P : (11.411 - [1.431) which can 
be ca st into the form of (11.491) : and assumes the concentration to remain uniform (so that 
(11.480 is suppressed, with constant in (I1.47P ). In the ma crosco p ic app roach to constitutive 
modelling, one instead arrives at equations of the form (11.471 - 11.490 by ansa tz, or by an 
exact description of a simplified system (such as the dumb-bell model of Section [1.3.50 . Fur- 
ther guidance comes from the generic constraints of translational and rotational invariance, 
and from Onsager reciprocity. The crucial advantage of the macroscopic approach is that 
it allows coupling between the flow and microstructural quantities such as 0, Qij and P{L) 
to be incorporated in a simple way. Another ingredient, almost alwa ys abs ent from micro- 
scopic models (though see Ref. [84]), is that operators such as G in (I1.49P should contain 
spatially nonlocal terms, which are needed to correctly describe the structure of spatially 
inhomogeneous {e.g., banded) flows [63]: 

DtW,, = VnVm, 0) + VV^W,, (1.51) 



1.5.1 Johnson-Segalman Models for Shear Thinning 

Within the framework just described, we now discuss some specific models of shear thinning. 
These are designed to reproduce, at the level of macroscopic modelli ng, th e nonmonotonic 
constitutive curve of the microscopic reptation-reaction model (Sec. II 4.3p . for which ho- 
mogeneous flow is unstable with respect to the formation of shear bands. From now on, 
we reserve the term "constitutive curve" for the underlying nonmonotonic relation between 
stress and strain rate, and use "steady state flow curve" for the actual stress/strain-rate 
relation measured in an experiment. In the banding regime, the former is unstable and gives 
way to the latter: the two only coincide when the flow remains homogeneous. 

The most widely used model was originally devised by Johnson and Segalman [85], and 
later extended by Olmsted et al [86] to include the spatially nonlocal terms needed to describe 
the struct ure o f the interface between the bands. For ce ba lance and incompressibility are 
given by ( 11.450 and ( 11.460 . The viscoelastic stress of ( 11.47P is assumed to depend linearly 
on the molecular deformation tensor Wij, and on the concentration via a modulus G: 



(1.52) 
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The deformation tensor Wij obeys diffusive Johnson-Segalman (dJS) dynamics as follows: 

DtW,, = a{DuWij + WuDi^) + {Wuni, - naWi^) + 2D,, - + VV^W,, (1.53) 

in which Dij and Qij are respectively the symmetric and antisymmetric parts of the velocity 
gradient tensor Kij = VjVi. (Note that the deformation tensor Wij used here differs from 
that in Section [1.3 51 by a trivial isotropic contribution.) For the moment we assume that 
is uniform so that G(0) and r((/)) are constants in space and time for any sample. 

S etting a = 1 and P = in (11.531) we recover Oldroyd B dynamics, as derived in Sec- 
tion 11.3.51 by considering an ensemble of relaxing dumb-bells undergoing affine deformation 
under shear. Although, as described there, Oldroyd B is th e most natural extension to 
nonlinear flows of the linear Maxwell model of Section ll.3.2[ its trivial constitutive curve 
(^xyii) = G'jT + r/^°^7 fails to capture the dramatic shear thinning and related flow instabil- 
ities seen in wormlike micelles. To allow for shear thinning, the Johnson-Segalman model 
invokes a 'slip parameter' a, with \a\ < 1. When \a\ < 1, the dumb-bells no longer deform 
affinely, but slip relative to the flow field. The resulting constitutive curve is then 

The viscoelastic contribution (first term) now shear-thins dramatically and is nonmonotonic, 
increasing as G'jr for small 7 before decreasing towards zero at higher shear rates (the 

maximum is at 'yry/l — a? = 1). In contrast, the Newtonian solvent stress always increases 
with 7. The overall shape of the constitutive curve thus depends on the relative strength of 
these contributions. For r]^°^ > Gt/8 the solvent dominates the viscoelastic stress to restore 
monotonicity (bottom dashed curve in Fig. 11.121) . For lower 77^°^ the negative slope survives 
over some range of shear rates (upper dashed curves in Fig. 11.12]) . Assuming that the modulus 
increases with concentration as G ~ 0^'^ and the relaxation time as r(0) ~ 0^'^ [32] , one 
obtains the full family of dashed curves (Jxyii, 0) of Fig. 11.121 [87]. As discussed in Sec. 11.5.31 
below, the spatial gradient terms in (11.531) are needed to confer a finite interfacial width 
/ oc VVr between the bands [63] . 

In work related to the above, Ref. [88], Cook and Rossi considered the dynamics of an 
ensemble of bead-spring dumb-bells of finite extensibility, subject to spring, Brownian and 
hydrodynamic forces, allowing for slip. They thereby derived a coupled set of equations of 
motion for the stress and the number density of dumb-bells. When the finite extensibility 
parameter tends to zero, these reduce to the Johnson-Segalman model. They studied the 
predictions of this model, including shear banding, in planar shear for a variety of boundary 
conditions. In Ref. [89], this work was generalised to Taylor-Couette fiow. 



Concentration Coupling 



Above we assumed the micellar concentration to remain spatially uniform. Generically, 
however, one expects concentration fiuctuations to be important in sheared multi-component 
solutions when different species have widely separated relaxation times [90-93], as seen 
experimentally in Refs. [94-97]. This was first explained by Helfand and Fredrickson [98] in 
the context of polymer solutions, as follows. Under shear, parts of a stretched polymer chain 
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Figure 1.12: Dashed lines: constitutive curves for the diffusive Johnson- Segalman model 
with G ~ 0^'^, r ~ 0^'^ and various 0. Circles: limits of linear stability in the dJS model. 
Squares: corresponding limits for the full dJS0 model with an experimentally realistic mi- 
cellar diffusion coefficient. (Diamonds and triangles are for an artificially reduced T>.) Inset: 
zoom on large 7. Figure reprinted with permission from Ref. [87]. 



(or micelle) in regions of low viscosity will, on relaxing to equilibrium, move further than 
parts mired in regions of high viscosity and concentration. A relaxing chain thus on average 
moves towards the high concentration region. This provides a positive feedback whereby 
chains migrate up their own concentration gradient, leading to fiow-enhanced concentration 
fiuctuations. 

In a remarkable paper, Schmitt et al [99] outlined the implications of this feedback mecha- 
nism for the shear banding transition. They predicted an enhanced tendency to form bands, 
together with the existence of a concentration difference between the bands in steady state. 
They futher noted that this difference would lead to a slight upward ramp in the stress 
"plateau" of the steady state flow curve. Subsequently, strongly enhanced concentration 
fluctuations were seen in the early-time kinetics of shear band formation [100] . Observations 
of a ramping stress plateau are now widespread; for example, see Ref. [67 . In Ref. [87], 
therefore, one of us proposed an extension to the dJS model, by combining the constitutive 
equation fll.531) with a two-fluid model for flow-concentration coupling. 

The basic assumption of this two-fluid approach [90-93] is a separate force-balance for 
the micelles (velocity w™, volume fraction 0) and solvent (velocity vf, volume fraction 1 — 0) 
in any fluid element. Any relative velocity = — vf (implicitly assumed zero in the 
ordinary dJS model) can then give rise to concentration fluctuations. The forces and stresses 
acting on the micelles are assumed as follows: (i) the usual viscoelastic stress G{(f))Wij; (ii) 
an osmotic force (j)ViSF/S(j) derived from a free energy F, leading to conventional micellar 
diffusion; (iii) a drag force C{4>)vl^^ impeding motion relative to the solvent with a drag 
coefficient (; (iv) an additional Newtonian stress 2(j)ri^^D^^ due to fast micellar relaxation 
processes such as Rouse modes, where Df-^ is the symmetric traceless part of the micellar 
strain rate tensor; and (v) a hydrostatic pressure. The solvent experiences (I) the usual 
Newtonian viscous stress; (II) a drag force (equal and opposite to that on the micelles) and 
(III) a hydrostatic pressure. Now adding (i)-(v), one obtains a force balance equation for 
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the micelles. Adding (I)-(III) we get the corresponding equation for the solvent. The sum 
of these two equations gives the force balance equation for the fluid element as a whole 



6F{<P) 



= V,G'(0)W^,, - 0V,^^^ + 2V,0r7P°i Af + 2V, (1 - 0)r7^°^ A° - (1.55) 

which replaces (11.451) . (We have redefined r]^°^ slightly here by explicitly extracting the 

prefactor 1 — 0.) Another combination of the two equations yields the relative velocity vl^^, 
which governs the concentration fluctuations: 



C(0) 

The evolution of the molecular deformation tensor Wij is prescribed by dJS dynamics, as 
before, but now with the velocity Vi in (11.531) re-interpreted as the micellar velocity vf^. 

The constitutive curves of this "dJS0 model" , which by definition describe purely homo- 
geneous flow, are the same as for the original dJS model. (Dashed lines in Fig. 11.121 ) The 
relevance of the new coupling is that any heterogeneity in the flow variables now affects the 
concentration field, and vice versa. As we will show below, this enhances the tendency to 
form shear bands (Sec. [T75^ and leads to a concentration difference between the bands in 
steady state (Sec. 11.5.31) . as earlier predicted in Ref. [99]. In the limit of large drag, ( ^ oo, 
concentration fluctuations are suppressed, and we recover the original dJS model. A two- 
fluid model of shear banding was independently developed by Yuan and Jupp in Ref. [101]. 
Here it was claimed that gradient terms in the free-energy functional for the concentration 
field are sufficient to give unique stress selection (see Sec. ll.5.3] below) for the banded state, 
without the need for gradient terms in the viscoelastic constitutive equation. (This is at odds 
with Ref. [68]; the discrepancy deserves further investigation, but might be attributable to 
the particular numerical scheme employed in [101].) 



A Simplified Scalar Model 



The dJS model is the simplest tensorial model to capture the negatively sloped constitutive 
curve of the full reptation-reaction picture. An even simpler model [102] neglects normal 
stresses and considers only the shear stress a = a^y and the shear component 7 = VyV^ of 
the strain rate tensor. It further chooses units that equate both r and G to unity, so that 
(jPy^ = Wxy The force balance and constitutive equations are then 

a = aP°l + r^«°l7 (1.56) 

r9taP°l = -aP°l + gijr) + fV^a^"^ (1.57) 
in which the choice g{x) = x/(l + x^) is made; this recovers a JS-like constitutive curve. 



1.5.2 Shear Banding Instability 



The above models each have a non-monotonic constitutive curve 17(7). The negatively sloping 
part of this is well known to be mechanically unstable [103]. Under conditions of imposed 
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shear rate, the system can recover steady flow by separating in to she ar bands, one on each of 
the stable, positively sloping branches (discussed fully in Sec. 11.5.31 below). In this section, 
we discuss the initial onset of banding in the unstable regime, by means of a linear stability 
analysis. 

For simplicity, we start by analyzing the simple scalarised model of fll.56[ [L57j) . Consider 
a homogeneous initial state at an applied shear rate 70. The viscoelastic shear stress is 0"^°^ = 
gi'joT) and the total shear stress a = gi'joT') + f]^°^o- Consider now a small heterogeneous 
perturbation away from this state, decomposed into Fourier modes with wavevectors confined 
for simplicity to the flow gradient direction y: 



t) = a^' + E cos (^) eMu^nt) (1.58) 



with a similar expression for 7. We assume flow between parallel plates aX y = and y = A 
with plate conditions of zero-slip for the velocity and zero-flux for the stress, 9^0"^°^ = 0. 
The Un are growth rates to be determined. If Un < OV_n, the flow is stable; Un > signifles 
instability and the onset of banding. Linearising fll.56p and (11.571) in these perturbations, 
and recognising that the total stress must remain uniform by force balance, we obtain 

= 5(tP°i + r7^°i57„ (1.59) 

c^r<5ar^ = -5crr^ + g'iioT)5% - PklSa^' (1-60) 

in which the wavevector is = rnr/A, and the function g{x) and parameter / are as deflned 
in (ll.57p . Combining these, 

oOnT = -^5'(7or) - I'kl (1.61) 

in which S' = g'ij^r) + 77^°^ is the slope of the underlying constitutive curve. When this is 
negative we have positive growth rates ci;„ > for wavevectors /c„ less than a large cutoff 
(scaling as 1/0- This signifles linear instability and the onset of shear banding. A related 
analysis, for the Doi-Edwards model for unbreakable polymers, can be found in Ref. [104]. 

The assumption above is that the system starts on the underlying homogeneous constitu- 
tive curve. In practice, it is impossible to prepare an initially homogeneous state within the 
unstable region. Instead, one performs a shear startup experiment at the desired shear rate 
7o. The linearisation must therefore now be done about this time- dependent startup flow, 
complicating the analysis somewhat [105]. However, the basic stability properties turn out 
to be unaltered from the simplifled calculation presented above. 

In the scalarised model, there is only one dynamical variable, crP°'. What about the dJS 
model of (11.531) ? In principle, tensoriahty now confers three dynamical stress components, 
dependent on three components of the strain rate, VjWj, with ij = xx,xy,yy. (We have 

suppressed concentration and set r(0) = 1 and G{(j)) = 1, so af°^ = Wij.) However, if we 
allow spatial variation only in the flow-gradient direction y, incompressibility conflnes the 
velocity to the flow direction x, and only one component VyVx = 7 of the strain-rate tensor 

is relevant. It further turns out that, of two linear combinations Y, Z of aP°i and aP°i, only 
Z is relevant: Y decouples and decays. The governing equations then reduce to the form 

^^y = < + (1-62) 
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Figure 1.13: Left: homogeneous flow. Riglit: shear banded flow. 

= fii, <, Z) + fy'alf (1.63) 
Td,Z = g{^,alf,Z) + l'V^Z (1.64) 

For stabihty purposes the terms can now be neglected since they act only to cut off 
any instability at high wavevectors [kl ^ 1) as seen above. Linearising these equations in 
small perturbations about a state of homogeneous flow on the underlying constitutive curve, 
we then obtain a quadratic equation for the growth rate: u;^ + hun + c = 0. For this to have 
at least one unstable root, Un > 0, we need 6 < or c < 0. For this model it can be shown 
that c = cS', where S' is the slope of the constitutive curve, and b and c are positive at all 
shear rates. Therefore, the condition for instability is S" < as before: homogeneous flow is 
unstable in the intermediate range of shear rates where the constitutive curve has negative 
slope (circles in Fig. 11.121) . 



1.5.3 Steady Shear Bands 



For an applied shear rate in the regime of decreasing stress, we have seen ab ove t hat homo- 
gene ous f low is unstable. For a constitutive curve of the shape shown in Fig. 11.121 (or indeed 
Fig. ll.lOl above for the reptation-reaction model, or Fig Jl. 141 below for dJS), this instability 
triggers formation of two bands of shear rates 71 and 72, one on each of the stable branches, 
with a flat in terface between bands whose normal is oriented in the flow-gradient direction 
y (Fig. I1.13p . The relative volume fractions (/, 1 — /) of the bands arrange themselves to 
match the spatially averaged shear rate 7 imposed on the cell as a whole. (It is this averaged 
quantity, namely the ratio of the velocity difference between plates to the width of the gap, 
that now appears on the abscissa of the experimental flow curve, (1(7).) As 7 increases, the 
width / of the high shear band increases at the expense of the low shear band. Force balance 
demands that the shear stress a^y is common to both the bands. Assuming that the nature 
of each band does not vary as their amounts change {i.e., neglecting concentration coupling), 
this gives a plateau in the observed flow curve (7(7). 

The scenario just described, and flrst proposed for micelles in [60], was conflrmed by 
explicit numerical calculation within the dJS model [106]. The resulting steady state flow 
curve is shown in Fig. 11.141 and indeed comprises two homogeneous branches (7 < 71 and 
7 > 72) connected by a plateau across the banding regime 71 < 7 < 72- A typical flow 
proflle in the banding regime is shown in Fig. 11.151 Note the smooth variation across the 
interface, which has a width / oc VVt set by the spatially non local (stress diffusion) term 
in (11.531) . Without this term, the interface would be an sharp discontinuity whose position 
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Figure 1.14: Constitutive curve (dashed line) and steady state flow curve (solid where differ- 
ent) in the dJS model (a = 0.3, r]^°^ = 0.05, G = r = 1). Figure reprinted with permission 
from Ref. [106]. 



is not reproducible between different startup runs [86]. The nonlocal term confers a smooth 
interface and a robust, reproducibly selected value o"sei of the plateau stress. The selected 
stress is then the only one at which a stationary front can exist between the bands [63]. 
We note that an alternative nonlocal JS model incorporating a higher order gradient of 
the deformation-rate tensor was offered by Yuan in [146] and solved numerically in two 
dimensions (the fiow/fiow-gradient plane). The model was shown to give a uniquely selected 
stress, although the model was not fully resolved in the sense that the width of the interface 
between the bands depended on the fineness of the numerical mesh. 

In recent years, the experimental evidence for shear banding in wormlike micelles has be- 
come overwhelming [34,55-58,67]. Reports of kinks, plateaus and non-monotonicity in the 
flow curve are now widespread, while spatially resolved NMR [55-57] and birefringence [58] 
data provide direct evidence for banding in both shear rate and microstructure. As noted 
previously, in some cases, the associated stress plateau is not perfectly fiat, but ramps up- 
wards from left to right. In a cylindrical Couette geometry, there will always be a small slope 
caused by a slight stress gradient (absent in the planar case of Fig. I1.13P causing the high 
shear band always to reside next to the inner cylinder. As this band expands outwards with 
increasing applied shear rate into regions of lower stress, the overall torque must increase to 
ensure that the interface between the bands stays at the selected stress asei- An alternative 
explanation of the upward slope, independent of cell geometry, is coupling between flow and 
concentration, as first discussed in Ref. [99] (see Sec. II. 5. 11 above). If a concentration differ- 
ence exists between the bands, the properties of each band must change as the applied shear 
rate is tracked through the coexistence regime, because material is redistributed between 
them as the high-shear band expands to fill the gap. This was confirmed by one of us in 
Ref. [68] by a numerical study of the dJS0 m odel. The stress now slopes upwards from left to 
right through the shear banding regime (Fig. ll.l6ll . The concentration and shear rate in each 
phase now define a family of tie li nes, one for each banded state, giving the nonequilibrium 
phase diagram shown in Fig. 11.171 

Beyond the Johnson- Segalman model, shear banding has also been studied in the Doi 
model of shear thinning rigid rods [107-109]. In this case, the relevant microstructural vari- 
ables in fll.47p are the nematic order parameter Qij and the concentration of rods 0. While 
this approach obviously ignores any effects of micellar flexibility, it takes a flrst step to incor- 
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Figure 1.15: Shear banded profile predicted by evolving the dJS equations in one spatial 
dim ension , y. The imposed shear rate 7 = 2.0, towards the left of the stress plateau in 
Fig. 11.141 G = T = 1. Figure reprinted with permission from Ref. [106]. 



porating orientational ordering, ignored by dJS0 and likely to be important in concentrated 
micellar solutions close to an underlying isotropic-nematic transition. Depending on param- 
eter va lues, the constitutive curve for homogeneous flow can now adopt either of the shapes 
in Fig. 11.181 In both cases, the two stable branches correspond to a low-shear isotropic band 
(branch I) and a flow- induced paranematic phase (branch II). For concentrations inside the 
zero-shear biphasic regime, branch II touches down to the origin to form the zero-shear ne- 
matic phase. Coupling between concentration and flow arises because more strongly aligned 
rods in the high shear band can pack more closely together, giving a higher concentration. 
In contrast, the Helfand-Fr edrick son coupling in the dJS0 model gives a less concentrated 
high shear band; recall Fig. I1.17[ 

The shape of the lower constitutive curve in Fig. 11.181 opens up a new possibility, shown 
by the vertical line: that shear bands can coexist at a common shear rate with a different 
value of the stress in each band. In a Couette cell, this corresponds to bands stacked in 
the vorticity direction (Fig. 11.191 right) and is called "vorticity banding". In contrast, the 
dJS model supports only common-stress banding (horizontal line in the (7, axy) plane), with 
the normal to the banding interface in the flow gradient directi on; th is is called "gradient 
banding". The latte r give s concentric bands in Couette flow (Fig. II 19l left), reducing to the 
arrangement of Fig. 11.131 in the limit of planar shear. 

In this section, we have explored steady shear-banded states within models whose under- 
lying constitutive curve comp rise s two stable branches separated by an unstable region of 
negative slope. In Sees. II. 5. 61 and II. 5. 71 below, we discuss the effects of (i) higher dimension- 
ality and (ii) more exotic constitutive curves. As we will see, either can give rise to unsteady, 
even chaotic, shear bands. First, however, we introduce some models of shear thickening. 
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Figure 1.16: Steady state flow curves in the dJS0 model for different values of the average 
concentration 0. Concentration coupling now confers an upward slope in the banding regime. 
Figure reprinted with permission from Ref. [68]. 
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Figure 1.17: Nonequilibrium phase diagram: tie lines show the two-phase (banding) regime. 
Different tie lines correspond to different values of the stress. Dashed line shows the spinodal 
limit of linear stability of homogeneous flow. Figure reprinted with permission from Ref. [68]. 
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Figure 1.18: Stress-strain-rate curves for the Doi model with different excluded volume 
parameters. The dashed line segments are unstable. The straight lines indicate possible 
coexistence between branches I and II under conditions of common stress (horizontal lines) 
or strain rate (vertical line). Figure reprinted with permission from Ref. [109]. 




Figure 1.19: Geometries for phase separation at common stress (left) or strain-rate (right) 
in a Couette rheometer. At a common stress (left) phases I and II have different strain rates, 
while at a common strain rate (right) they have different stresses. Here z is the vorticity 
axis, X is the flow direction, and y is the flow gradient axis. Figure reprinted with permission 
from Ref. [109]. 
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Figure 1.20: Shear thickening constitutive curve. Figure reprinted with permission from 
Ref. [80]. 




Figure 1.21: Bare (hght hne) and final (heavy) constitutive curve in the CHA model. The 
region of instability a~ < a < is shown. Figure adapted from Ref. [110]. 
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1.5.4 Simple Models for Shear Thickening 

As discussed in Sec. 11.4.61 above, a window of shear thickening is seen in many micellar 
solutions, for volume fractions around the onset of viscoelasticity [12,74-78]. After a period of 
prolonged shearing, an initially inviscid fluid undergoes a transition to a notably more viscous 
state or even a long-lived gel, with shear banding often implicated in its formation [77]. 

There is no consensus on the microscopic origin of this phenomenon (see Sec. I1.4.6p . 
Despite this, several macroscopic features can be captured within a phenomenological ap- 
proach [80,81] that couples flow to a generalised gelation transition by allowing a mixture of 
two species A (sol; concentration 0^ = 1 — 0) and B (gel; concentration 0^ = 0) to slowly 
inter-convert under the influence of shear: 

at0 = i?(0,7)+WV (1-65) 

i?(0,7) = |7|(l- 0)02-^0 (1.66) 

In the absence of shear, = 0_b relaxes to zero, leaving pure A. In this way, B is identifled 
as the shear-induced phase (gel). The model of [80] ignores normal stresses and considers an 
additive shear stress a = aA + ctb, with both contributions presumed Newtonian 

a = [(1 - 4>)r]A + (Pvb] 7 (1-67) 

Hence, all other variables are deemed fast in comparison to the structural variable 0, an 
assumption that only strictly becomes valid in the limit of vanishingly slow inter-conversion. 
Depending on the ratio f]B/VA = c, the model can actually capture either shear thinning or 
shear thickening. Her e we f ocus on thickening, c > 1, for which the underlying constitutive 
curve is shown in Fig. 11.201 

As discussed in Ref. [80], for an applied shear rate 7 = F in the non-monotonic regime, 
the system can in principle choose between homogeneous states on branches 1 or 3 (circles in 
Fig. I1.20p or it can gradient-band between these branches at a selected stress a* (horizontal 
line). But when the model equations are evolved numerically at imposed shear rate, gradient 
banding is not seen. The system instead always chooses homogeneous flow: on branch 1 
below a critical shear rate 7 = 7* (Fig. 11.201) . and on branch 3 for 7 > 7*. The sharp vertical 
jump between branches 1 and 3 at 7* leaves a range of stresses that is unattainable under 
conditions of controlled shear rate. This range is instead accessed by controlling the stress, 
and marks a regime of vorticity banding. 

A different model of shear thickening was devised by Gates, Head and Ajdari (CHA) in 
Ref. [110]. This has an unspecifled, slowly evolving structural variable, but now with explicit 
(though still relatively fast) dynamics for the stress. The distinguishing feature of this model 
is that, while the "instantaneous" constitutive curve following fast stress equilibration at flxed 
structure is non-monotonic (Fig. ll.2Tl) . the slow structural evolution restores monotonicity in 
the true long term constitutive curve. (In contrast, in Fig. I1.20l the non-monotonicity results 
directly from the structural evolution.) The short-term tendency to form shear bands is thus 
frustrated by the long-term stru ctural evolution, leading to oscillatory or chaotic dynamics. 
We return to this model in Sec. 11.5.61 below. 

It appears that the microscopic mechanisms involved in shear banding models differ from 
case to case, and often remain poorly understood. Nonetheless, the macroscopic phenomenol- 
ogy has many universal features, including kinks, plateaus or non-monotonicity in the flow 
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curve. In fact, there are just four fundamental types of flow curve in shear banding sys- 
tems: thinning vs thickening, and gradient vs vorticity banding. Ref. [Ill] collects these 
curves, discusses their dependence on concentration, and explains how a non-equilibrium 
phase diagram can be reconstructed from a family of such curves. 



1.5.5 Other Related Models 



Here we discuss some other approaches to modeling the rheology of wormlike micelles. In an 
article this length we have no room to discuss all of the many nonlinear rheological models 
that were developed without micellar systems specifically in mind. Of these, alongside the 
JS-type models described extensively above, we mention only the Giesekus model [112] of 
shear thinning polymeric fluids, based on the concept of a deformation dependent tensorial 
mobility. This was first applied to micelles in Ref. [44] ; it shows a plateau in the constitutive 
curve, but lacks the non-monotonicity required to give a true banding instability. Any ability 
to fit the measured flow curve for shear-banding micelles is therefore somewhat fortuitous. 

Among micelle-inspired approaches, Manero and co-workers [4,113,114] developed a sim- 
ple model coupling the evolution of the stress to that of the underlying fluid structure. For 
simplicity, they represented the structure by a single "fluidity" parameter = l/rj, assumed 
to be the reciprocal of the fluid viscosity. (Here is distinct from the micellar concentra- 
tion, for which the symbol was used above.) According to Ref. [114], possible microscopic 
interpretations of l/<y9 include the number of bonds, links or entanglements in a network. No 
distinction is made between the viscoelastic and solvent stresses: instead, the total stress is 
directly assumed to obey an upper convected Maxwell constitutive equation 

+ 7^iKu(Tlj + (TiiKji) = ^ (1.68) 

in which is the symmetrised shear rate tensor and Gq is the plateau modulus. The 
fluidity is assumed to evolve as 

^ = fl_f. + k[y]{(poc - v)(^ijDji (1.69) 

in which yjQ and (poo are the steady state fluidities in the limits of zero and inflnite shear 
rates; A is a structural relaxation time for the build up of structure in zero shear, while /c[7 
is a function allowing for the breakdown of structure under shear. For ip^ > ipQ, the mode 
captures shear thinning. (Predictions for shear thickening, yjoo < V^o? are not discussed here.) 

In the simplest case k is assumed constant, independent of shear rate. Here the con- 
stitutive curve is always monotonic, ranging from Newtonian for ip^o = (po through simple 
shear thinning for cp^o > </3o > to yield stress behaviour for ip^o > (po = 0. This version of 
the model therefore lacks a shear banding instability. Ref. [114] discusses its predictions for 
various rheological tests in homogeneous shear, including step strain, step stress, shear-rate 
jumps, and stress sweeps. In particular, signiflcant hysteresis is predicted in upward followed 
by downward stress sweeps. To allow for shear banding, in Ref. [113] the rate of structure 
breakdown is assigned an additional dependence on shear rate by taking 

fc[7] =fco(l + /ii|7|) (1.70) 

For large enough values of /xi and 'Poo/'^o, the homogeneous constitutive curve_js non- 
monotonic, assuming the basic form shown for the Johnson- Segalman model in Fig. II. 141 The 
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model thus captures a shear banding instabihty. The predictions of Eqns. fll.68[ 11.691 11.701) 
for the ste ady sh ear banded state are discussed in Ref. [115]. In this case, nonlocal terms 
(recall Sec. 11.5.31) are not invoked to capture a uniquely selected banding stress. Instead, the 
authors appeal to extended irreversible thermodynamics to posit the existence of a Gibbs 
potential that can be used to calculate state selection. However, no deeper justification 
was offered for this assumption. Shear rate profiles for banded states were not detailed in 
Ref. [115], but in the absence of nonlocal terms, these presum ably a re discontinuous across 
the band interface, in contrast to the models discussed in Sec. 11.5.31 

Goveas and Pine [116] developed a simple model of shear-thickening in a wormlike micellar 
solution. Above a critical stress, the solution is assumed to undergo a reaction that produces 
an insoluble gel phase. Competing against this is the destruction of gel by, for example, 
peeling at the interface between gel and solution. The balance of these determines the 
relative volume fractions of gel and sol, and therefore the position of the interface between 
the two. The following assumptions were made: (i) a constant viscosity for the Newtonian 
sol phase, (ii) no flow in the gel phase, (iii) a constant rate of gel destruction, and (iv) a rate 
of gel production proportional to the product of the stress and the surfactant concentration 
in the sol phase. Given these, the model predicts a re-entrant region in the flow curve 
that is only accessible under controlled stress conditions. For controlled shear rate, the 
measured flow curve is predicted to be discontinuous. These results agree qualitatively with 
the experiments of Pine on shear thickening micelles [77]. 

In Ref. [67], Porte, Berret and Harden argued for a description of shear-banding involving a 
modified thermodynamic transition as distinct from a mechanical instability. On this basis, 
they developed a thermodynamically inspired model based on a potential energy Fs(7s), 
where 7^ is an internal structural variable related to the local recoverable elastic strain. The 
local stress is assumed to be directly set by this strain as 

dFs 

with the following dynamics assumed 

at a^fs 

in which /i is a mobility. Here the first term on the right is a dissipative term; the second is 
a driving (reactive) term. Solving Eqn. 11.721 in homogeneous steady shear gives 

a = -^7 (1-73) 

in which l//i(7s) is the non- Newtonian viscosity. As shown in Ref. [67], the homogeneous 
constitutive curve (7(7) is non-monotonic, leading to banding, whenever the relation cr(7s) is 
itself non-monotonic. Because the underlying nonlinearity stems from the dissipative term 
—^{'-ys)dFs/d'ys, this approach effectively considers banding as a flow-induced parameter shift 
for an underlying thermodynamic transition, rather than the result of an intrinsically rheo- 
logical instability. While this model may indeed be relevant to micellar systems close to an a 
true thermodynamic transition (for instance to a state of nematic order) [67], the suggestion 
in Ref. [67] that a robust stress plateau and sigmoidal kinetics [64] are incompatible with a 
purely mechanical instability was disproved by subsequent work on the dJS model [86, 117]. 



More generally, although a distinction is frequently made between 'thermodynamic' and 
'mechanical' instabilities, it is clearly possible to construct models including ingredients of 
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both. This was done in Ref. [68], by couphng the JS model to a nearby demixing instabil- 
ity. Here, the matrix that determines the linear stability of homogeneous flow contains a 
concentration subspace (which can be subject to a zero-shear demixing transition); and a 
'mechanical' subspace (which is subject to a banding instability even without concentration 
coupling). Furthermore, the coupling between these subspaces mixes the two instabilities, 
and a smooth crossover can be achieved between the two by varying the model parameters. In 
the same way, the distinction between 'mechanical' banding instabilities and shear-induced 
isotropic- nematic instabilities is likely to be less clear cut than was once thought. 

Dhont introduced a phenomenological model of shear banding in Ref. [118], as follows. 
At finite Reynolds number in planar shear, with spatial variations only in the flow-gradient 
direction y, the Navier- Stokes equations is 

dvjy^t) ^ dcr{y,t) 
^ dt dy ^ ^ ^ 

where v is the fluid velocity in the flow direction and p is the mass density. The total shear 
stress cr is assume to respond adiabatically, such that 

^{y, t) = viiiv, miy, t) - «:(7(y, t))^^^ (1.75) 

The flrst term on the right describes the nonlinear homogeneous constitutive relation between 
stress and shear rate. The second term encodes the spatial gradients needed to give unique 
state selection in the banding regime. Within this term, k, called the "shear curvature 
viscosity" in Ref. [118], is further assumed to obey lim^^oo '^(t) = 0. In steady state, the 
total shear stress must be a constant, independent of y (see Eqn. I1.74p . Denoting this by 
(Jstatj Dhont showed that the banded state obeys 

d7 [^^h(7) - ^stat]/«:(7) = (1.76) 

where Oh^y) is the part of the total stress obtained by discarding the gradient term in 
Eqn. fll.75p . If k were constant, fll.76p would reduce to the equal area (Maxwell) construction 
that determines equilibrium phase coexistence. Note, however, that this is somewhat special 
property of a scalar model in which normal stresses are neglected. (In a tensorial model, 
even with constant k, no equal area construction would be possible in general.) 

Dhont then performed a linear stability analysis for fluctuations ab out a homogeneous 
initial state of shear rate 70, to study the onset of banding, as in Sec. 11.5.21 above. For a 
fluctuation of wavelength kn, he found a growth rate 



UJr. 



^70 



+ «:(7o)fcr, 



^ (1.77) 
P 



As in Eqn. f ll.6ip above, instability is predicted when the slope of the homogeneous consti- 
tutive curve is negative. In contrast to fll.6ip . however, it predicts a maximum growth rate 
at a non-zero wavevector 



max 



rffTh(7o)/ci7o 



due to a competition between diffusion of (conser ved) m omentum, governed by fll.74p . and 
the gradient terms in (11.750 . This contrasts with fll.6ip . where the competition is between 
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the rate at which the (non-conserved) stress evolves, and the gradient terms. For reahstic 
micellar parameters and rheometer gap sizes, the latter appears to us more appropriate. 

In a separate study, Dhont and Briels [84] developed an expression for the stress tensor 
in inhomogeneous suspensions of rigid rods, in terms of the flow velocity and the probability 
density function for the position and orientation of a rod. By explicitly allowing for large 
spatial gradients in the shear rate, concentration, and orientational order parameter, this 
approach could potentially be applied to shear-banding systems. 



1.5.6 Temporal Instability 



The studies discussed so far capture the basic tendency of wormlike micelles to undergo a 
transition to shear banded flow. However most of them fail to address recent reports that the 
constitutive response to steady mechanical driving is intrinsically unsteady in some regimes. 
In such cases, the stress response to a constant applied strain rate (or vice versa) does not 
settle to a constant value. Instead it shows sustained periodic oscillations [119-123] or erratic 
behaviour suggestive of low-dimensional chaos [57, 124- 126]. Such long-time unsteadiness is 
distinct from the early-time instability discussed in Sec. 11.5.21 the latter merely provides the 
initial trigger for banding, of some sort, to occur. 

Hydrodynamic instabilities (Taylor-Couette, Rayleigh-Benard, turbulence, etc.) have 
long been studied in simple liquids [127], where they stem from the nonlinear inertial 
term (/WjVjfj) in the Navier-Stokes equation. For the flows of interest to us, however, 
the Reynolds number is virtually zero, rendering this term negligible. The observed com- 
plexity must therefore instead arise in the constitutive nonlinearity inherent to the rheology 
of the micelles themselves, e.g., through the coupling between microstructure and flow. The 
term "rheo logical chaos" [110] (or "rheochaos") has been coined to describe this behaviour. 
Irregular signals have also been reported in other complex fluids, include so-called 'onion' 
surfactant phases [128, 129] and concentrated colloids [130]. In many cases, evolution of 
the microstructure in concert with the rheological signal has been explicitly observed via 
birefringence imaging [120,131], light scattering [128], or spatially resolved NMR [126]. 

A crucial question is whether these instabilities are spatiotemporal or purely temporal in 
character. In wormlike micelles they most often arise close to the banding regime, suggesting 
the spatiotemporal evolution of a heterogeneous {e.g., banded) state. Indeed, early optical 
experiments on wormlike micelles showed a temporally oscillating state comprising spatially 
alternating turbid and clear bands [120,121]. More recent advances allowing spatially and 
temporally resolved measurements of velocity profiles have unambiguously revealed fluctu- 
ating shear bands in both wormlike micelles [126] and multilamellar onion phases [129,132]. 
However, spatial observations have been made in just a few of these cases, so that in others, 
the question remains open. For example, a nearby banding instability could feasibly play 
a role in triggering temporal rheochaos, but with banding itself narrowly averted such that 
the system stays homogeneous, as discussed in Ref. [110]. For systems close to the nematic 
transition, another possibility is the purely temporal director chaos captured theoretically in 
models of nemato-dynamics [133,134]. In this section, we therefore discuss purely temporal 
instability i n spat ially homogeneous models, before proceeding to the full spatio-temporal 
case in Sec. I1.5.7[ We address thinning and thickening systems in turn. 
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Figure 1.22: Dynamical phases of steady and transient states at the I-N coexistence temper- 
ature. The sohd hne is the border between the in-plane orbits tumbling (T), wagging (W) 
and aligning (A); the dashed and dotted lines delimit the regions where the out-of-plane 
orbits kayaking-tumbling (KT) and kayaking- wagging (KW), respectively, exist. Here 7, 
Aeq and Ak denote dimcnsionless shear rate and tumbling parameters of the Erickson-Leslie 
theory of nematic hydrodynamics, with Ak = AeqOeq. Figure reprinted with permission from 
Ref. [140]. 




Figure 1.23: Trajectories of the eigenvector corresponding to the largest eigenvalue of UiUj, 
after the transient. Orbits are plotted over the unit plane, (a) - (d) correspond to increasing 
values of the shear rate. Figure reprinted with permission from Ref. [133]. 
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Shear Thinning: Rigid Rods 

Models of rigid rods have been widely studied in the context of liquid crystalline polymers. 
They capture an isotropic to nematic (I-N) transition, and ma ke pre dictions for director 
dynamics in the nematic phase under shear. As noted in Sec. 11.5.31 above, such models 
obviously ignore any effects of micellar flexibility or scission; but they do take a first step 
to incorporating orientational ordering, relevant for some micellar systems that are close to 
the nematic transition. The interplay of shear banding with the I-N transition was studied 
in detail (for rigid rods) in Refs. [109]; a recent review can be found in Ref. [135]. 

The studies in question, following Hess [136] and Doi [107], consider a population of rods, 
with orientation vectors Ut chosen from a distribution ip{ui,t). This is assumed spatially 
homogeneous. Taking account of macroscopic flow, excluded volume effects and thermal 
agitation, the evolution of ip is specified via a Fokker-Plank equation. To solve this numer- 
ically, one must first project it onto a finite number of degrees of freedom. One method 
is to expand ip in a truncated set of spherical harmonics, giving a set of coupled ordinary 
differential equations for the time-dependent expansion coefficients [137]. An alternative is 
to construct the second order orientation tensor Qij oc (uiUj)^, and project the dynamics 
onto it via a closure approximation [138] to get an evolution equation for Qij{t). 

The resulting numerical predictions for director dynamics in the nematic phase under 
shear can be briefly summarised as follows. Studies that confine the director to lie in the 
fiow/flow-gradient plane predict a sequence of transitions from "tumbling" through "wag- 
ging" to "flow-aligning" with increasing shear rate [137], for suitable values of a "tumbling 
parameter" A. (For other values of A, flow-alignment occurs at all shear rates.) In the tum- 
bling and wagging regimes the director executes periodic motion in the flow/flow-gradient 
plane. Studies generalised to allow out-of-plane director components [139] predict a richer 
dynamics, including new periodic regimes of "kayaking-tumbling", "kayaking- wagging" , as 
well as chaos characterised by a positive Lyapunov exponent and a fractal correlation di- 
mension [133,140]. Both intermittency and period-doubling routes to chaos are seen. The 
various regimes, calculated via Qij [140] are summarised in the phase diagram of Fig. 11.221 
Director trajectories, calculated using spherical harmonics [133], are shown in Fig. 11.231 



Shear Thickening: CHA Model 

As noted in Section [1.5.41 above. the CHA model of [110] couples the viscoelastic shear stress 
to a slowly evolving structural variable. The dynamics is defined as follows: 

a = 7 - R[ai) - A(T2 (1-79) 

with the structural evolution modelled by "retarded stresses" 

ai{t)= f Mi{t-t')a{t')dt' for i = l,2 (1.80) 

The Mi{t) are memory kernels, each havin g an in tegral of unity. In the absence of relaxation, 
the first term on the right hand side of (11.79^ causes the stress to increase linearly with 
straining (a Hookean solid with a spring constant of unity). The second and third terms 
respectively capture nonlinear and linear stress relaxation. 
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Figure 1.24: Time series of the stress at 7 = 18.49, 7 = 30, 7 = 33.38 (bottom to top) in 
the CHA model. CHA parameters: A = 20, T2 = 10 and R{a) = O.Qa^ - 3.3a^ + ba. Figure 
adapted from Ref. [110]. 

In the simplest version of the model. Mi is chosen to be a delta function such that 
(Ti{t) = a{t) is actually unretarded. For simplicity, the memory kernel M2 is chosen to 
be exponential, M2(t) = T2^exp (— t/r2). In this case, the system can be rewritten as two 
coupled differential equations in the stress a and structural variable m = a2'. 

d = 7 — i?(cr) — Am ; T2m = —{m — a) (1-81) 

In steady state at a given applied shear rate 7, we find the relation 

7 = i?(or) + Aa (1.82) 

When inverted, this defines the constitutive curve (7(7) of Fig. 11.21] Its two components -R(cr) 
and Act stem respectively from rapid nonlinear stress relaxation on a timescale t ~ = 
0(1) and retarded linear relaxation on a timescale t = T2 ^ 1. Thus -R(cr) represents an 
"instantaneous" constitutive relation, describing the relaxation of stress at fixed structure. 
The much slower structural relaxation eventually recovers the full curv e RM + Act. The 
interesting case arises when R{a) + Act is monotonic but R{a) is not (Fig. 11.21]) . The system 
then exhibits a shear banding instability, at short times, in the regime where R'{(t) < 0. If 
the linear contribution Act is sufficient retarded {t2 large), it fails to overcome this instability, 
despite the monotonic constitutive curve. Accordi ngly, the long-term dynamics of the model 
remain unsteady in a region a~ < a < (Fig- 11 ■241) . The dynamical system defined by 
(11.811) undergoes a Hopf bifurcation at ct+ and a~ (Fig. I1.2ip . signifying the onset of finite 
frequency sinusoidal oscillations with an amplitude varying as I7 — 7c|^''^- 

Chaos requires a ph ase space of dimensionality at least three; it cannot occur in the 
dynamical system (11.811) . Without invoking flow inhomogeneity (which gives infinite dimen- 
sionality), sufficient dimensions can be achieved by assuming cti = n to be retarded as well as 
CT2, with Mi(t) = rf^ exp (— t/ri). The underlying constitutive curve is again given by (ll.82p . 

In harmony with the simpler version of the model, one takes ri ~ 1 <^ r2 and consider the 
situation where monotonicity of the constitutive curve is restored only via the more retarded 
relaxation term. In the unstable regime, one now finds a period-doubling cascade leading to 
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temporal chaos. However, the physical interpretation of the second retarded stress is unclear 
and the model therefore remains strongly empirical in nature. 

The reader may be struck by the naive character of the CHA model when compared 
with the fully tensorial descriptions outlined in the preceding section for nematics liquid 
crystals. (Admittedly, the latter field has had an extra half-century to develop its equations!) 
Note however that the scalar treatment of shear stress is not quite as restrictive as it may 
appear. Specifically, it does not rule out the existence of normal stresses, but does assume 
that the shear stress has an autonomous dependence on strain rate in simple shear flows. 
This holds, for instance, in the upper convected Maxwell model despite the latter's fully 
tensorial character. The development of tensorially convincing equations for shear-thickening 
materials remains a topic for future study, but so far attention has focussed instead on 
relaxing the assumption of spatial homogeneity of the flow. This we address next. 



1.5.7 Spatiotemporal Instability; Rheochaos 

Above we have discussed the unsteady rheological response of models with purely temporal 
dynamics. Such models assume from the outset that the sample remains homogeneous, 
with each point in space following an identical evolution in time. As noted above, however, 
reports of unsteady dynamics in wormlike micelles are most common close to or inside the 
shear banding regime. In such spatio-temporal description is essential, to allow for an 

evolving state that is heterogeneous {e.g., banded) at any instant. Indeed, recent experiments 
have unambiguously revealed fluctuating shear bands in both wormlike micelles [126] and 
multilamellar onion phases [129,132]. We now turn to models of spatiotemporal dynamics, 
considering thinning and thickening systems in turn. 



Shear Thinning: A One-Dimensional Model 



As explained previously, homogeneous flow on the negatively sloping branch of t he con sti- 
tutive curve is unstable with respect to the formation of shear bands [103]. In Sec 11.5.2] we 
further saw that coupling between the flow and an auxiliary variable such as concentration 
can enhance this instability, causing it to extend into regions of positive constitutive slope 
(Fig. I1.12p . In Ref. [141], one of us exploited this fact to construct a model in which the 
high shear band is itself destabilised, leading to unsteady banding dynamics. 

The model is defined as follows. We work in one spatial dimension, the flow-gradient 
direction y, with a velocity Vi = v{y)6ix and shear rate 7(2/) = dyv{y). Normal stresses 
are neglected, and the total shear stress is assumed to comprise additive viscoelastic and 
Newtonian components. At zero Reynolds number, a{t) must be uniform across the gap: 

a(t) = aP°i(y,t)+rt(2/,t) (1-83) 
For the dynamics of the viscoelastic component, we use the scalar model of Sec. 11.5.11 

a^poi^.^ + sIilM + paj, (1.84) 

T[n) r(nj ^ 

with a relaxation time r that now depends on a structural variable n = n(y,t), according 
to T{n) = ro(n/?T,o)". As before, g{x) = x/{l + x'^) is chosen to ensure a region of negative 
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Figure 1.25: Intrinsic constitutive curves for differing degrees of coupling between flow and 
micellar length, a) Weak coupling, giving the standard coexistence of stable low and high 
shear bands (A and B); b) moderate coupling; c) strong coupling. Squares show Hopf 
bifurcations. The thin black lines delimit the periodic orbit of the local model at fixed a. 
Figure reprinted with permission from Ref. [141]. 
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Figure 1.26: Homogeneous dynamics for fixed stress, log{a) — —0.301, for different initial 
conditions. Figure reprinted with permission from Ref. [141]. 
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constitutive slope. The auxiliary variable n is taken to represent a nonconserved quantity. 
For definiteness it is identified as the mean micellar length (previously denoted by L) although 
there could be other candidates for its interpretation. Coupling of n to the flow is completed 
by assuming that it evolves with its own relaxation time r„, distinct from r(n): 



Here the coupling term N{x) = no/(l + x^) represents (say) shear-induced scission. 

As discussed in Sec. 11.5.21 above, a state of initially homogeneous flow on the underlying 
constitutive curve is unstable (growth rate u;„ > 0) if 6 < or c < in the dispersion 
relation + bun + c = 0, where c = cS' and S' is the slope of the constitutive curve. The 
coefficients b and c depend on the model parameters and the applied shear rate. As in the 
dJS model, it can be shown that c > 0; this model therefore shows the familiar banding 
instability when 5" < 0, at intermediate shear rates. It also shows a new instability, 6 < 0, 
not seen in dJS. This destabilises the high shear branch to a degree that depends on r„, 
eventually terminating in a Hopf bifurcation shown by squares in Fig. I1.25[ 

The resulting model's homogeneous dynamics are explored in Fig. 11.261 which shows 
parametric phase portraits 7(t), n{t) at a fixed value of the stress a for various initial 
conditions. Here, any possibility of spatial structuring has been artificially suppressed so 
that 7 and n can depend only on time. The new instability at high shear rates corresponds 
to an unstable focus in the 7, n plane, associated with a pair of complex conjugate eigenvalues. 
Trajectories originating near this point spiral outwards to attain a stable limit cycle. 

Turning to the full spatio-temporal dynamics at imposed shear rate 7 = / dy'^{y,t), for 
small r„, we find stable shear bands, connecting points A and B in the constitutive curve 
of Fig. 11.25b . For higher r„, the h igh-shear band is unstable, leading to unsteady dynamics 
of the banded state; see Fig. 11.271 where several regimes are evident. At low applied shear 
rate a thin pulse or band of locally high shear ricochets back and forth across the cell. (A 
thin fiuctuating high shear band, away from the rheometer wall, was seen experimentally 
in Ref. [57].) At larger shear rates, we find two or more such pulses. For two pulses (not 
shown), we typically find a periodically repeating state with the pulses alternately bouncing 
off each other (mid-cell) and the cell walls. Once three pulses are present {e.g., 7 = 7.0), 
periodicity gives way to chaotic behaviour. At still higher mean shear rate, 7 = 19.2, we 
find regular oscillations of spatially extended bands pinned at a stationary defect. The local 
shear rates span both the low and high shear constitutive branches. (Oscillating vorticity 
bands were seen experimentally in Ref. [131,142].) For the intermediate value 7 = 11.35 
we find intermittency between patterns resembling those for 7 = 7.0 and 7 = 19.2. Finally 
for 7 = 23.0, 31.0 we find oscillating bands separated by moving defects, with the fiow now 
governed only by the high shear constitutive branch: in each band, the shear rate cycles 
round the periodic orbit of the local model (Fig. 11.261) . 

For different t:„ we find other interesting phenomena [141]. For example, for weaker 
instability (r^ = 0.13) at low applied shear rates we see a high shear band that pulsates in 
width while adhering to the rheometer wall, or meanders about the cell (Fig. 11.281) . The 
former behaviour resembles that sometimes seen in micelles [56,143] and onion phases [144]. 
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Figure 1.27: Right panels: spacetime plots showing the time evolution of the shear rate 
profile for different average shear rates (as identified in the left panels). The data is for the 
ID model with parameters set at = 0.145, V = 0.0016. The space coordinate < y < 1 
runs left to right and time 120 < t < 140 runs bottom to top; the shading denotes local shear 
rate with regions of high shear rate showing as dark zones (as per the grey-scale bar at the 
bottom of each panel). Taking a horizontal slice through the spacetime plot at a given time 
shows as dark the regions of high shear rate; a slice higher up the plot shows these regions 
at a later time. Hence a localized band of high shear rate, moving across the system with 
velocity v, shows up as a grey stripe of slope v on the spacetime plots. Left panels: time 
series showing the corresponding stress vs. time. The average shear rate 7 is denoted 7 in 
the main text. Figures adapted from Ref. [141] and reprinted with permission. 
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Figure 1.28: Shear rate intensity plots for r„ = 0.13, D = 0.0016. Figure reprinted with 
permission from Ref. [141] 



Shear Thinning: Higher Dimensional Model 

In Sec. ll.S.ll we introduced the diffusive Johnson-Segalman (dJS) model and in Sec. I1.5.3l we 
discussed its predictions for steady shear banded states. Those calculations were restricted to 
ID variations only (in the flow gradient direction y) and implicitly assume that the interface 
remains perfectly fiat at all times. Recent experiments suggest more complex interfacial 
dynamics [126], and a natural question is whether these ID states remain stable in higher 
dimensions, or whether they give way to more complex behaviour. This question was recently 
addressed by one of us in a study comprising two separate stages: first, a linear stability 
analysis of an initially ID banded state [106], and then a numerical study of nonlinear 
interfacial dynamics [145]. We summarise each of these in turn. 

Linear stability analysis: As discussed in Sec. 11.5. H a ID flow-gradient calculation of 
planar shear within the dJS model predicts t he st eady state flow curve of Fig. I1.14[ with 
a corresponding banded profile shown in Fig. 11.151 To study the linear stability of this ID 
"base" profile with respect to small fluctuations with wavevector in the {xz) plane of the 
interface, we linearise the model for small perturbations (lower case) about the (upper case) 

base profile: ^{x,y, z,t) = $(y) + (f)q{y) exp^ujqt + iq^x + iqzz). The vector $ comprises 
all components $ = (WafiiVa)-, the pressure being eliminated by incompressibility. This 
gives a linear eigenvalue equation. We are interested in the eigenvalue ti;max(q) with the 
largest real part, 3f?co'max(q): in p articular, whether it is stable, SfJcUmax < 0, or unstable, 
J^i^max > 0. We find (Fig. I1.29p that at any q^, JJwmax increases with decreasing values 
of / = \JV>T^ which sets the width of the interface between the bands. For small enough 
I the dispersion relation is positive over a range of wavevectors, rendering the ID profile 
unstable. This applies to shear rates across most of the stress plateau of the flow curve, and 
is furthermore robust to variations in the JS parameter a. Because the / values accessed here 
(I r-u \ — lOyum for a 1mm rheometer gap) are even larger than those expected physically 
(/ ~ lOOnm) this study suggests that, experimentally, the entire stress plateau region will 
be unstable to perturbations away from a flat interface between shear bands. 

Nonlinear interfacial dynamics: Once the undulations attain a finite amplitude, nonlinear 
effects become important, and the linear calculation then breaks down: one must perform 
a full nonlinear study of the dJS m odel in the flow/flow-gradient (xy) plane [145]. The 
results are summarised in Fig. I1.30[ At large values of /, for which the dispersion relation 
of the linear analysis is negative at all wavevectors, the ID base profile remains stable as 
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Figure 1.29: Real part of the eigenvalue of the most unstable mode; a = 0.3, rj = 0.05, 
7 = 2.0, R eynolds number p/r] = 0. The data for / = 0.01 correspond to the base profile in 
Fig. I1.15[ Symbols: data. Solid lines: cubic splines. Figure reprinted with permission from 



expected: the interface stays flat at all times. For smaller values of /, the dispersion relation 
is positive over a window of qx (Fig. 11.291) : at fix ed /, the number N of linearly unstable 
modes increases with the system size A^,. (Fig. I1.30|) . For small // A^, ju st inside the unstable 
regime, the ultimate attractor comprises a travelling wave (Fig. II 3ip . The wall-averaged 
shear stress is constant in time, with a value W^y^ss that depends on / and A^. and is slightly 
higher than the selected value W^"^ of the ID calculation. For l/A^ values deeper inside 
the unstable regime in Fig. 11.301 we see a new regime in which the travelling wave now 
periodically "ripples" (Fig. 11.321) . The corresponding wall-averaged stress W^y is periodic in 
time, with variations of the order of one percent, and an average value larger than the ID 
selected stress W^'^. The interface height h{x, t) is shown as a white line in Fig. I1.32d . 

The results just discussed were obtained in numerical runs starting from an initial state 
comprising two adjacent shear bands separated by a flat interface. Shear startup from rest 
can give more complicated results, allowing multiple interfaces to form, with complex dynam- 
ical interactions suggesting low dimensional chaos. As discussed in Ref . [145] , cell curvature 
(almost always present experimentally) is likely to eliminate this degeneracy, restoring the 
single-interface scenario described above. 

Note that in earlier work. Yuan et al [146] used a mixed Lagrangian-Eulerian algorithm 
to evolve a related non local JS model in 2D planar shear. Contrary to Ref. [145], they 
apparently found a stable interface. However, the banded profiles reported in Ref. [146] were 
averaged along the flow-direction. In a ny case, rather short cell lengths A^^^ were used, likely 
to be in the stable regime of Fig. I1.30[ 
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Figure 1.30: Nonequilibrium phase diagram showing regimes of nonhnear interfacial dynam- 
ics in the dJS model together with the numbers of hnearly unstable modes. [L^ = in 
the present notation.) The slip parameter a = 0.3 and applied shear rate 7 = 2.0. Figure 
reprinted with permission from Ref. [145]. 



k « 

Figure 1.31: Greyscale of order parameters for travelling wave in the (x, y) plane for / = 0.015, 
Aj. = 6, and upper wall velocity V = 'jAy = 2 to the right. (Notation differs from the text: 
= Wij.) Figure reprinted with permission from Ref. [145]. 



Shear Thinning: Other Relevant Approaches 



Here we touch on two other theoretical approaches to spatio-temporal rheochaos in generic 
viscoelastic fluids, which may be relevant to micelles. The first generalises models of ne- 
matodynamics (Sec. 11.5.61) to allow spatial structuring, and finds spatio-temporal rheochaos 
setting in via intermittency [147,148]. Although the existing study neglects backfiow effects, 
a hydrodynamic equivalent could be relevant in concentrated wormlike micelles close to the 
nematic transition, where the highly sheared band may have a high degree of nematic or- 
dering. However, any direct link with wormlike micelles remains to be established; for more 
details see a recent review [135]. 

The second approach is based on the concept of "elastic turbulence" . For simple Newto- 
nian liquids, it has long been known that the nonlinear inertial term of the Navier Stokes 
equation can cause smooth, laminar flow to destablise at high Reynolds number, giving way 
to more complicated flow profiles [149,150]. For complex fluids, including wormlike micelles. 
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Figure 1.32: Rippling wave at / = 0.005, A^. = 4, 7 = 2. Greyscale of y). Upper wall 

moves to the right. White line in d): interface height. Figure reprinted with permission from 
Ref. [145]. 



inertia is usually negligible. However, analagous instabilities (which are somewhat distinct 
from those involving shear-banding as discussed above) can arise directly from constitu- 
tive nonlinearity. As reviewed in Ref. [151], this opens up the possibility of purely elastic 
instabilities that occur even in inertialess flows at zero Reynolds number. 

For example, the Oldroyd B model discussed in Sec. 11.3.51 predicts purely elastic insta- 
bilities triggered by hoop stresses in both curved Couette [152] and torsional [153] flows. 
Such phenomena have been observed experimentally in model "Boger" fluids, which com- 
prise high molecular weight polymers in viscous solvents, in both curved Couette [154] and 
plate-plate [155] flow. The basic observation is that, above a critical shear rate, the lam- 
inar base flow destablises, accompanied by significant increase in mean shear stress. In a 
remarkable development, a transition series linking the base flow at low shear rates to fully 
developed elastic turbulence [156-158] at high shear rates was uncovered. This involves mul- 
tiple dynamic states including stationary rings, competing spirals, multi-spiral chaos, spiral 
bursts and eventually fully developed elastic turbulence. Planar shear lacks hoop stresses, 
so the linear viscoelastic instabilities discussed above are absent. Nonetheless, a nonlinear 
equivalent is still predicted to arise, since fluctuations cause the originally parallel stream- 
lines to become curved, and so subject to instability at nonlinear order [159]. The relevance 
to micelles of these generic findings on elastic instabilities remains to be explored. 
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Figure 1.33: Nonequilibrium phase diagram of the shear-thickening model of Ref. [161] when 
Ts and (a) are varied. Black triangle: chaotic points, black circle: periodic point. Three main 
regimes are observed: (O) oscillating shear bands, (T) travelling shear bands, (C) chaotic 
regions. The outer dashed line is the linear stability limit. Figure reprinted with permission 
from Ref. [161]. 



Shear Thickening 

In Sec. 11.5.61 above, we discussed the temporal dynamics of the CHA model of shear thick- 
enin g, in w hich the instantaneous constitutive relation at fixed structure is nonmonotonic 
(Fig. 11.211) . This creates a short-term tendency to form bands of differing shear stress, coex- 
isting in the vorticity direction at a common value of the shear rate. However this short term 
tendency is opposed by the long term structural evolution. In a full spatio-temporal scenario, 
the interplay of these two effects can give rise to shear banding with complex dynamics. 

The CHA model of Sec. 11.5.61 was recently extended to allow for spatial structuring [160, 
161], allowing the stress to vary along the vorticity direction: a(t) — > a{z,t). A spatial 
gradient term PV^a was added, in line with t he dis cussion of Sec. ll.5.3l above. For simplicity, 
only (72 was taken as retarded, by extending (11.811) as follows: 

&{z,t) =^{t) - R{a) ~ Xm + VV'^a (1.86) 

m{z,t) = -^^^^ (1.87) 

As before the instantaneous nonlinear relaxation term was chosen as 

R{a) = ao - ba^ + ca'^ (1.88) 

giving in linear response the Maxwell time tm = l/o. Note that the shear rate 7(t) in (11.861) 
is now uniform: the moving wall of the rotor imposes the same velocity for all heights z, so 
that 7(z, t) = 7(t) only. The model is studied under conditions of imposed torque, i.e., under 
an imposed value of the spatial mean of the stress (a). The two main control parameters 
are taken to be (a) and the ratio Ts/tm- Depending on the values of these parameters, the 
model shows a rich variet y of o scillatory and chaotic banding dynamics, as summarised in 
the phase diagram of Fig. 11.331 
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Figure 1.34: Various types of spatiotemporal rheochaos observed in the shear thickening 



model of Ref. [160, 161 
spacetime plots (as exp 
with permission. 



Left: time series for a local shear rate and for stress. Right: 
ained in Fig. 11.271) . Figure adapted from Ref. [161] and reprinted 
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In the regime w here t he structural evolution is much slower than the stress relaxation 
(marked 'O' in Fig. I1.33p . the model predicts oscillating shear bands. Varying the imposed 
stress (a) along any horizontal line of fixed Ts/tm in this regime, the associated waveforms 
range from simple to very complex. Near the middle of the line one sees simple 'fiip-fiopping' 
bands, with the cell divided equally between a high-stress and a low-stress band, the identities 
of which repeatedly switch with a period of order the structural time Tg. Moving slightly 
off-centre along the line, the interface between the bands now adopts a zig-zagging motion, 
superposed on the fiip-fiopping motion just described. Near the edges of the lin e, ve ry 
complex oscillations are seen. The regime 10 < Ts/tm < 10^ marked 'T' in Fig. 11.331 is 
characterised by a periodic nucleation of shear bands, which then cross the cell with roughly 
constant velocity. In the two disconnected pockets marked 'C, at moderately long values of 
Ts relative to tm, and strongly off-centered values of (a), the complex oscillations of regimes 
'O' and 'T' give way to true rheochaos, characterised b y a p ositive Lyapunov exponent. 
Various examples of rheochaotic flows are shown in Fig. 11.341 This model also admits an 
interesting low-mode truncation in which only the lowest two non-homogeneous Fourier 
modes are retained for each of 7 a nd n, giving 4 modes in all [161]. The basic structure 
of the phase diagram of Fig. 11.331 was found to be preserved by this truncation, showing 
that rheochaos is robust within the model, and is not dependent on the presence of sharp 
interfaces between the bands. Within this truncation, the chaos in the pockets 'C was found 
to set in via a classical period doubling scenario. 



1.5.8 Rheochaos: Relation of Macroscopic Theories to Experiment 

From the preceding sections it is clear that a wide range of oscillatory, irregular and chaotic 
behaviour can be found within relatively simple models of both shear-thinning and shear- 
thickening fluids, once structural memory is allowed for. Some of these models were directly 
intended to address giant micelles; others were not. 

Over the past several years, there have been a range of experimental studies, many of them 
referred to above, in which broadly similar flow behaviour was reported in micellar systems. 
There is an emerging consensus that in micelles, spatiotemporal rather than purely temporal 
instability is the norm. However, it is in the very nature of unstable and chaotic systems that 
detailed prediction of particular trajectories is a near-futile task; these trajectories depend 
on both model parameters and initial conditions in an erratic fashion. The more important 
task, in confronting theory and experiment, is to develop robust tests of whether the physical 
content of the models is supported by or explains the experimental data. Such tests are likely 
to concern scenarios (such as whether the route to chaos is by period doubling, or some type 
of intermittency) rather than trajectories. At the time of writing, such detailed comparative 
work has barely begun. 

In many related flelds, including work on inertial hydrodynamic instabilities in Newtonian 
fluids, the task is hampered by the difficulty of acquiring clean experimental data that can 
expose the route to chaos and/or the low-dimensional stucture of the underlying attractors on 
flrst entering the chaotic regime. Interestingly, careful recent work shows such measurements 
to be possible with very high precision in certain micellar systems; see for instance [126,162]. 
The recent data of Ganapathy and Sood [163] deserves special mention. It shows a transition 
to chaos via a type II intermittency route, and gives direct evidence for coupling between 
concentrations and shear stress. This raises the hope that an increasing body of high quality 
experimental data will soon allow discrimination between the various theoretical models 
outlined above, and lead to better ones being formulated. 
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1.6 Summary and Outlook 



In this article, wc have reviewed theoretical modelling efforts that address the rheology of 
giant micelles. In the well-entangled regime, the linear viscoelastic spectra, often close to 
pure-Maxwell in character, are robustly explained by the reptation-reaction model which 
couples micellar kinetics to the tube dynamics of entangled objects. The occurrence of a 
single relaxation time, despite the presence of exponential polydispersity in micellar lengths, 
arises from the rapidness of micellar reactions on the time scale of stress relaxation and the 
consequent averaging over micellar reaction dynamics of the local rates for stress decay. 

In addressing the concentration dependence and ionic strength dependence of the Maxwell 
time, the reptation-reaction model has had relatively mixed success; however this is in part 
attributable to difficulties connecting the input parameters of the model (micellar end-cap 
energy, reaction rates, etc.) to quantities more directly measurable in experiment. Setting 
aside these difficulties, once one allows for micellar branching — which, perhaps surprisingly, 
assists rather than impedes stress relaxation — the model does seem to account for the main 
observational trends in the regime of linear viscoelasticity. 

For well-entangled micellar materials, shear-banding is routinely observed in the nonlin- 
ear flow regime. This was an early, robust prediction of the reptation-reaction model, which 
gives a nonmonotonic dependence of shear stress on strain rate in steady shear flows. (The 
model can also rationahze the ever-increasing normal stress observed as the mean shear.) 
The original reptation-reaction model was however based on a 'ffist generation' description 
of the tube which also predicts nonmonotonic flow curves in monodisperse, unbreakable 
polymers. In the past decade, 'second generation' tube models, incorporating the concept 
of convective constraint release have been developed. So long as the phenomenological pa- 
rameter describing the strength of this effect is neither too large nor too small, these 'second 
generation' models can restore the observed monotonicity of the flow curve for unbreakable 
polymers, while retaining nonmonotonicity for micelles. The primary rheological predictions 
of the reptation-reaction model thereby remain intact. 

Beyond explaining the basic tendency to form shear bands, it has not yet proved practi- 
cal to relate the rather complex nature and dynamics of such bands to micellar constitutive 
models directly. One must instead turn to macroscopic models which couple features of the 
reptation-reaction approach to spatially varying order parameters. These models currently 
offer several competing descriptions of a wide range of phenomena involving unsteady or 
chaotic banded flow, many of which were observed in recent experiments. Further detailed 
experimental work is highly desirable so as to enable better discrimination between the vari- 
ous models on offer. In particular one would like to clarify the interdependence between the 
following effects: local stress relaxation; concentration coupling (which causes the micellar 
volume fraction to vary between bands); stress diffusion (which can control the physics of 
band interfaces); normal stress effects at such interfaces; curvature of streamlines (causing 
'purely elastic' instabihties) ; and structural memory. The last term refers to the evolution 
of internal degrees of freedom, such as the mean micellar length, on a time scale that can be 
at least as large as the characteristic time for stress relaxation. More than one of this list 
of physical ingredients could be implicated in the destabilization of steady shear bands, but 
we are not sure yet what is the balance of power among them. Additionally, the flows seen 
in entangled micelles can sometimes include bands of extremely high shear rate, where there 
could be very strong molecular alignment at the micellar scale. This could allow molecular 
physics of a nonuniversal kind (beyond that captured by coupling to liquid crystalline order 
parameters and/or shear-rate dependent micellar reaction rates) to influence the macroscopic 
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rheology. This could occur in regimes where the average shear rate would not lead one to 
suspect any such effect. 

The position is even more complicated for weakly entangled micellar systems in the con- 
centration window close to the onset of viscoelasticity. Here, drastic shear-thickening can be 
seen, such that at high flow rates the system becomes much more strongly entangled (or at 
least, much more viscoelastic) than at rest. The resulting shear-induced gel-like phase can be 
long lived, suggesting a role for structural memory in this case also. Further evidence comes 
from the fact that the thickening behaviour is often accompanied by shear banding and/or 
unsteady flows, with complex sample-history dependences in some cases. Microscopic mod- 
elling of this esoteric but fascinating regime remains speculative, and offers limited guidance 
to formulating macroscopic models. However, one candidate, which can plausibly explain 
various of these experimental features, is a scenario in which the presence of micellar rings 
comes to dominate the rheology. Whether or not this eventually proves to be the correct 
picture, models in which the local Maxwell time is coupled to slowly evolving structural 
variables (such as the mean micellar chain length, or a gel fraction parameter) may offer a 
promising way forward in addressing the bafflingly complex experimental physics of giant 
micellar systems close to the onset of viscoelasticity. 

On the experimental side also, the onset regime has tended to be neglected, perhaps be- 
cause the fully entangled case is more relevant for applications. (The latter, not reviewed 
here, include personal care products, oilbore fluids, and creation of new biocompatible ma- 
terials [4].) Another long-neglected area has been the study of entangled micelles in strongly 
elongational flows. Study of such flows has, in recent years, been a major driver in the de- 
velopment of improved constitutive models for unbreakable polymers, where these flows are 
crucial in process areas such as fibre-spinning. Clarifying the similarities and differences be- 
tween micelles and conventional polymers in such flows could help to guide future modelling 
efforts for both classes materials. 
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